JLAB-THY-07-651, MIT-CTP-3841, TUM-T39-07-09 



Nucleon Generalized Parton Distributions from Full Lattice QCD 

Ph. Hagler 

Institut fiir Theoretische Physik T39, Physik-Department der TU Miinchen, 
James-Franck-Strafie, D-85747 Garching, Germany 

W. Schroers* 

John von Neumann-Institut fiir Computing NIC/DESY, D-15738 Zeuthen, Germany 

J. Bratt, J.W. Negele, and A.V. Pochinsky 
Center for Theoretical Physics, Massachusetts Institute of Technology, Cambridge, MA 02139 

R.G. Edwards and D.G. Richards 
Thomas Jefferson National Accelerator Facility, Newport News, VA 23606 

M. Engelhardt 

Physics Department, New Mexico State University, Las Cruces, NM 88003-8001 

G.T. Fleming 

Sloane Physics Laboratory, Yale University, New Haven, CT 06520 

B. Musch 

Institut fiir Theoretische Physik T39, Physik-Department der TU Miinchen, 
James-Franck-Strafie, D-85747 Garching, Germany 

K. Orginos 

Department of Physics, College of William and Mary, Williamsburg VA 23187-8795 
and Thomas Jefferson National Accelerator Facility, Newport News, VA 23606 

D.B. Renncr^ 

Department of Physics, University of Arizona, 1118 E 4th Street, Tucson, AZ 85721 

(LHPC Collaboration) 
(Dated: December 9, 2009) 

We present a comprehensive study of the lowest moments of nucleon generalized parton distri- 
butions in Nf = 2+1 lattice QCD using domain wall valence quarks and improved staggered sea 
quarks. Our investigation includes helicity dependent and independent generalized parton distribu- 
tions for pion masses as low as 350 MeV and volumes as large as (3.5 fm)'', for a lattice spacing of 
0.124 fm. We use perturbative renormalization at one-loop level with an improvement based on the 
non-perturbative renormalization factor for the axial vector current, and only connected diagrams 
are included in the isosinglet channel. 
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I. INTRODUCTION 

Generalized parton distributions (GPDs) P, [3, H, Q play a vital role in our understanding of the structure of the 
nucleon in terms of the fundamental building blocks of QCD, the quarks and gluons. Before the advent of GPDs, 
fundamental questions as to the origin of the spin of the nucleon, the decomposition of the nucleon total momentum, 
and the distribution and density of the nucleon constituents in position and momentum space seemed to be largely 
unrelated. In some cases, it was even unclear how to formulate these questions in a theoretically sound way and 
how to measure the underlying observables experimentally. With the introduction of GPDs, it is possible not only to 
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FIG. 1: GPDs as part of a scattering amplitude. 

give precise definitions to quantities, such as the quark and gluon angular momentum contributions to the nucleon 
spin and the probability densities of quarks and gluons in impact parameter space @, but also to unify and 
extend the successful concepts of parton distribution functions (PDFs) and form factors. Nucleon generalized parton 
distributions are experimentally accessible in deeply virtual Compton scattering of virtual photons off a nucleon and 
a range of other related processes 0, H, ■ Since these reactions involve in general convolutions of GPDs over the 

longitudinal momentum fraction a;, which makes it difficult if not impossible to map them over the whole parameter 
space, the most stringent quantitative information on GPDs currently comes from quark PDFs and nucleon form 
factors (Tl| . 

Complementary to experimental efforts, lattice QCD offers a unique opportunity to calculate x-moments of GPDs 
from first principles. The first investigations of GPDs including studies of the quark angular momentum contributions 
to the nucleon spin have been presented by the QCDSF collaboration in quenched QCD T3| and by LHPC/SESAM 
in Nf = 2 lattice QCD [l^. Lattice results on nucleon GPDs published since then have provided important insights 
into the transverse structure of unpolarized nucleons [3 , the lowest moments of polarized and tensor GPDs , 
and transverse spin densities of quarks in the nucleon [l7l . [Tsj l . With the exception of several initial studies [l^, HO] , all 
previously published lattice results on GPDs have been obtained from calculations in a two-flavor " heavy pion world" 
with pion masses in the range of 550 to over 1000 MeV. In this work, we improve on previous studies by presenting a 
comprehensive analysis of the lowest three moments of unpolarized and polarized GPDs in Nf = 2 + 1 lattice QCD 
with pion masses as low as 350 MeV and volumes as large as (3.5 fm)'^. 

The paper is organized as follows. We begin with an introduction to the calculation of moments of GPDs in lattice 
QCD in section II. Section III describes the hybrid approach of using domain wall valence quarks with 2+1 flavors 
of improved staggered sea quarks. In section IV we present our numerical results for the generalized form factors, 
including a discussion and interpretation of the quark orbital angular momentum and the transverse nucleon structure. 
Chiral extrapolations of selected lattice results to the physical pion mass are presented in section V. Conclusions are 
given in the final section. 



II. LATTICE CALCULATION OF MOMENTS OF GENERALIZED PARTON DISTRIBUTIONS 

Generalized parton distributions determine off-forward matrix elements of gauge invariant light cone operators 
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where x is the momentum fraction, n is a light cone vector and T ~ fi or T — fi^^. The twist-2 tensor GPDs [21 1 
related to F = rif^a^^ , j — \,2 are not studied in this work. The four independent twist-2 unpolarized and polarized 
generalized parton distributions H, E, H and E arc defined by the parametrizations 

n A 

{P\A'\0^{x)\P,A) ^ mH{x,U) + ^{{^'yn)E{x,U) , (2) 

and 

~ 7? • A ~ 

(P',A'|a/^,(:E) |P,A) = {{MH{x,U) + -^{{75))E{x,U), (3) 

2m 

where we use the short-hand notation ((F)) = U{P' , A')rU{P, A) for products of Dirac spinors U, and where A — 
P' — P, t — and £, = —n ■ A/2. In Eqs. ([2]) and ([3]) we suppress the dependence of the GPDs on the resolution 
scale fj? . An illustration of the GPDs parametrizing the lower part of the handbag diagram is given in Fig. ([1]). The 
momentum fractions x and ^ both have support in the interval [— 1,+1]. Depending on x, there are three kinematic 
regions, which offer different interpretations for the GPDs. For x G 1] and x G [—1,—^], the GPDs describe the 
emission and reabsorption of quarks and anti-quarks, respectively. In the case that x lies in the interval [— ^,^], they 
describe the emission of a gg-pair. 
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In our lattice calculations, we do not work directly with the bi- local operators in Eq. ([T]) but instead consider 
moments, defined by the integral J^^ dx x"'-^ f{x), of the operators in Eqs. (j2l3p . leading to towers of symmetrized, 
traceless local operators 

where [75] denotes the possible inclusion of the corresponding matrix, the curly brackets represent a symmetrization 

over the indices fj,i and subtraction of traces, and D~ ^/'2,{D — D)- We relate nucleon matrix elements of the tower of 
local operators in Eq. ([4]) to a;-moments of the twist-2 GPDs. To this end, we parametrize off-forward matrix elements 
(P', A'|C''i'^i- -''"^|P, A) in terms of the generalized form factors Ani{t), Bni{t), C„o(*), and Bni{t). Apart from 

potential difficulties related to lattice operator mixing for higher moments n, lattice measurements of the operators 
in Eq. ^ become increasingly noisy as the number of derivatives increases, and we therefore restrict our calculations 
to n < 3. The parametrization of nucleon matrix elements of Eq. ^ in terms of generalized form factors (GFFs) for 
n = 1, 2 and 3 reads 

Zra m 
2m 

+ A^''iA^^((7^-^>»A32(t) + :r^A^^^A''^((a^^}"))A„i?32W, (5) 

zm 

for the unpolarized case, and 

{P'\0^l \P) = ((7^^75Mio(i) + ^A^^ ((75))Sio(i) , 

(p'|0{mM.}|p) = p{Mi((^P2}^^))22o(t) + :^A{^^P^^>((75))B2o(t), 

(p'|0{MiP2M3}|p^ = P^^^P^M(7^^>75M30W + ^A^^'P^-P^''H{^5))B3o{t) 

'° 2m 

+ A{^iA^^((7^^>75))l32(i) + r^A{'^iA'^-A'^3}((^g))P32W, (6) 

2m 

for the polarized case. Here and in the following we set P = {P' + P)/2. 
Using Eqs. (|2l3l5l6p it is easy to show that Mellin-moments of the GPDs, 



77"(^,i) = j dxx'^-^H{x,^,t), £;"(^,t)EE J dxx'^-^E{x,^,t), 

^"(^,i)= / dxx"-^H{x,^,t), E"{C,t)= f dxx"-^E{x,^,t), (7) 
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are given by polynomials in the longitudinal momentum transfer ^ and the GFFs. For the lowest three moments, the 
corresponding relations read 
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(8) 

The aim of our calculation is to extract the GFFs as functions of the momentum transfer squared, t, from nucleon 
two- and three-point-functions as described below. Once the GFFs have been obtained, the complete i^-dependence 
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FIG. 2: Connected (a) and disconnected (b) diagrams in unquenched lattice QCD with an operator insertion at r and finite 
momentum transfer A. 

of the moments of the GPDs is directly given by Eqs. ([8]). Let us note that the Mellin-moments in Eq. (O are taken 
with respect to the entire interval from x = —1 to +1. Following our discussion below Eq. ([3]), we find that the 
moments of the GPDs at ^ = correspond to sums and differences of contributions from quarks q and anti-quarks q. 
For example, 



dxx''-\Eq[x,Q,t) + {-iy'Eq{x,Q,t)) , 
dxx''-\Hg(x,0,t) + (-l)("-i)i7^(x,0,i)) 
dxx''-\Eg{x,0,t) + (-l)("^i)£;^(a;,0,t)) . 



(9) 



Such a simple decomposition is not possible for non-zero longitudinal momentum transfer ^ 7^ 0. We denote the forward 
limit values of the moments of and H in Eq. ® by (a;"-i), = H^{0, 0) = A«o(0) and (a;"-i)A9 = ^^,"(0, 0) = l^o(O), 



where 



and 



)a(/ correspond to the moments of unpolarized and polarized quark parton distributions. 



Below, we give a brief summary of the methods and techniques used to extract moments of generalized parton 
distributions in lattice QCD. For details, we refer the reader to [3, Hsj. As usual, the matrix elements, Eqs. (|5l6p . 
are calculated from the ratio of nucleon three-point and two-point functions: 



C^P' (r, P) = J2 (runpoi),,. {Nk (r, P)N, (r,,, , P)) 



C^p'(t, P', P) = (rpoi),fe (iV,(r,„k, P')OiT, A)iV,(r,,e, P)) 



(10) 



where Funpoi = (1 + 74)/4 and Fpoi = (1 + 74)(1 + i7573)/2. The nucleon source, N{t, P), and sink, N{t, P), create 
and annihilate states with the quantum numbers of the nucleon. To maximize the overlap with the ground state, we 
used the smeared sources given in [2^. The three-point-function C^^{t,P',P) with the operator insertion at r is 
illustrated in Fig. ^ in terms of quark propagators, showing examples of connected and disconnected contributions 
in an unquenched lattice calculation. 

Using the transfer matrix formalism, we can rewrite Eqs. pop to obtain 



C2P*(r,P) = e-^°(^»"-"-) {Z{P)Z{P)f'^ -I- higher states. 



E^{P) 



(11) 



higher states , 



(12) 
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where the factors a,b, . . . represent the prefactors (including Dirac-matrices) of the corresponding GFFs A{t), B{t), . . . 
in the parametrizations in Eqs. (|5l6p . transformed to Euchdean space. Higher states with energies Ei > Eq in Eqs. (fTT|) 
and 112) are suppressed when Tgnk — t :s> l/(£'i — Eq) and r — Tgrc ^ ^/{Ei — Eq). 

In order to cancel the exponentials and Z-factors in Eq. (|12p for zero and non-zero momentum transfer A, we 
construct the ratio of two- and three-point-functions 

C^P^^snk - r + T,,c, P) C^P*(t, P') C^P^T^nk, P') ] . . 

C2pHrsnk-r-fT,,e,P')C2pHr,P)C2pt(r3„k,P)J ■ ^ ^ 

For an operator-insertion sufficiently far away from the source and the sink in the Euclidean time direction, higher 
states are negligible, and the ratio Ro{t, P' , P) exhibits a plateau in t. We finally average over the plateau region 
from Tniin to Tmax to obtain an averaged ratio Rq^P' , P). On a finite periodic lattice with spatial extent Lg, three- 
momenta are given by P = 2t: / {aLa)n with integer components rii — —Ls/2, .... La/2, and for the nucleon energy we 

use the continuum dispersion relation P4 = \/ w? + P^. Therefore, the discrete lattice momenta result in a finite set 
of values for the momentum transfer squared t which can be realized in our calculation. 

In order to obtain symmetric and traceless operators O'^'^^^^ "^, we have to choose specific linear combinations 
of the indices. For the diagonal operators typical examples are O^^^ = {O^^ -\- — O^^ — 0'^'^)/2 and 0^^2 — 
^Qi22 0133 „ 20'^'^^) /V2, where 05*=^ belongs to the 3-dimensional irreducible i?(4)-representation ^ for n = 2 
and 0"^2 is ^ member of the 8-dimensional representation r}^'' for ri = 3 [IJ] ■ The sets of operators Of we are using 
are the same as in [l^. Altogether, there are 9 linearly independent index combinations for n = 2 and 12 for n = 3. In 
order to be able to compare our results with experiment, the operators have to be renormalized and transformed to the 
MS-scheme at a renormalization scale n^. In general, operators mix under renormalization, and the renormalization 

matrix is non-diagonal. We will denote the renormalized operators in the MS-scheme by 0"'^^ — Z^Oj'. Some 
details concerning the renormalization procedure and numerical results for the renormalization constants will be 
discussed at the end of the next section. 

Based on the renormalized operators, we compute the averaged ratio Rq{P' , P) and equate it with the continuum 
parametrization in terms of the GFFs given in Eq. (jl2p . This is done simultaneously for all momentum combinations 
P and P' corresponding to the same momentum transfer squared t and all contributing symmetric and traceless 
operators Qf-^^^ giving a finite set of linear equations 

RoAPiPi) = ciiA{t) + ci2B[t) + . . . , 

RoAP2,P2) = C2lA{t) + C22B{t) + . . . , 
Po,™(P3,P3) = C3iA(t) + C32PW + . . . , 

(14) 

where (Pj — Pj)^ = t for all j = 1, 2, 3 . . .. The coefficients c,j in Eqs. (fT4|) only depend on the nucleon mass m and 
the momenta P, P' and are calculated from the traces in Eq. p^ . Finally, the set of equations ([H]), which in general 
is overdetermined, is solved numerically to extract the GFFs. The statistical errors for the GFFs are obtained from 
a jackknife analysis. 



RoiT,P',P)^ 



C^J\t,P',P) 

C2Pt(T.„k,P') 



III. LATTICE SIMULATION USING DOMAIN WALL VALENCE QUARKS WITH STAGGERED SEA 

QUARKS 

Since calculations at physical quark masses are prohibitively expensive with current algorithms and machines, 
we have used dynamical quark configurations at the lightest masses available, and where feasible, have used chiral 
perturbation theory to extrapolate to the physical mass. Staggered sea quarks with the Asqtad improved action were 
chosen because the computational economy of staggered quarks enabled the MILC collaboration to generate large 
samples of configurations at low masses on large spatial volumes [2^, [2^, which they freely made available to the 
lattice community. 

Chiral symmetry is crucial for avoiding some operator mixing, convenient for operator renormalization, and valuable 
for chiral extrapolation. Furthermore, the four tastes associated with staggered fermions immensely complicate 
calculating operator matrix elements in nucleon states. Hence, we chose a hybrid action utilizing chirally symmetric 
valence domain wall fermions (DWF) on an improved staggered fermion sea. Although this hybrid scheme breaks 
unitarity at finite lattice spacing, given the arguments that the valence and sea actions separately approach the physical 
continuum limit [27j . we expect that the hybrid action also approaches the physical continuum limit. Furthermore, 
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partially quenched mixed action chiral perturbation theory calculations are now becoming available for quantitative 
control of the continuum limit. We also note that hybrid actions have been successfully used in other contexts where, 
for example, the NRQCD action for valence quarks was combined with improved staggered sea auarks(28l| and was 
successful in predicting mass splitting in heavy quark systems. 

In our calculation, we used MILC configurations [l^ both from the NERSC archive and provided directly by the 
collaboration. We then applied HYP-smearing [sO] and bisected the lattice in the time direction. We have chosen 
gauge fields separated by 6 trajectories. Furthermore, we alternate between the first temporal half (time slices to 
31) and the second temporal half (time slices 32 to 63) on successive gauge configurations. In these samples we did 
not find residual autocorrelations. The scale is set by the lattice spacing a — 0.1241 fm determined from heavy quark 
spectroscopy [3l| with an uncertainty of 2%. 

Domain- wall fermions [H, [H, introduce an additional fifth dimension, L5 . They preserve the Ward-Takahashi 
identity [35l | even at finite lattice spacing in the limit L5 ^ cx). At finite values of L5 a residual explicit breaking 
of chiral symmetry is still present which can be parameterized by an additional mass term in the Ward-Takahashi 
identity (Sa . [37l |. In our calculations, we have kept this additional mass, (am)rcs, at least an order of magnitude 
smaller than the quark mass, (am)^^^ [l9j. To the extent that (am)res is negligible, perturbative renormalization 
of C[-yr,] is independent of 75 in the chiral limit and the non-perturbative renormalization of quark bilinear currents 
yields the same renormalization coefficients for the axial and the vector currents in the chiral limit j38j . 

We now consider the parameters entering the DWF action. The domain wall action realizes chiral symmetry by 
producing right-handed states on one domain wall that decay exponentially away from the wall and exponentially 
decaying left handed states on the other wall. To the extent that no low eigenmodes associated with dislocations (or 
rough fields) destroy the exponential decay and that the fifth dimension L5 is large enough, chiral symmetry will be 
nearly exact and (aTO)ros will be small. HYP smearing was essential to reduce the effect of low eigenmodes, but it 
is still necessary to use spectral flow to determine a value of the domain wall mass, M, for which the density of low 
eigenmodes is as small as possible. This was done on an ensemble of test configurations, with the result that we use 
M = 1.7. As discussed below, L5 was then tuned to keep (am)res below 10% of the quark mass and to have negligible 
effect on our lattice observables. Finally, the quark mass was tuned to produce a pion mass equal to the Goldstone 
pion mass for the corresponding MILC configurations. 



A. Tuning the fifth dimension 

The extent of the fifth dimension, L5, has been adjusted such that the residual mass, {am)-[cs is at least an order 
of magnitude smaller than the quark mass itself. This tuning is most relevant at the lightest quark mass since in 
that case the computational cost is largest and thus our tuning should be optimal. In addition, the breaking of chiral 
symmetry is also expected to be the largest and the resulting L5 provides a minimum value needed for our calculations 
at the higher masses. The residual explicit chiral symmetry breaking characterized by (0771)1-08 is obtained from f36'| 

^>^Al^2m,J^^{x) + 2J^x), (15) 

where 

J^^{x)^m,,,J^{x), (16) 

which holds up to 0{a?). 

We have run simulations using two samples of 25 configurations with volume Q. — 20"^ x 32: three degenerate 
dynamical Asqtad quarks with bare masses (a777,)Asqtad,sea _ q q5q (denoted as "heavy" and corresponding to m^r ~ 760 

MeV) and two plus one quarks with masses (am)^'"!''''*'"'''' = 0.010,0.050 (termed "light" and corresponding to 
TO^ ~ 350 MeV). The corresponding bare DWF masses have been adjusted to (oto)^^^ ~ 0.0810 and 0.0138 for the 
heavy and hght cases, respectively, cf. Sec. lIIIBl 

The resulting residual masses obtained from Eqs. and are plotted in Fig. ([3]). In the light quark case, 
— 16 just fulfills our requirement, while in the heavy quark case L5 — 16 more than satisfies it. This confirms our 
expectation that the value of L5 chosen at the lightest quark mass sets the lower limit for the other masses as well. 

One quantitative check that L5 = 16 is adequate is provided by the dependence of masses on L5 as shown in 
Figs. ([1]), ([5]), © and ([7]). The leading effect of m^cs is to shift the quark mass, so that when L5 is sufficiently large 
that this is the only effect, rn^ (x {niq + rrires)- Figs. ([4]) and ([5]) show the difference in the ratio m'^/{mq + mj-cs) 
at a general value of L5 and at L5 = 16, and indicate that the difference is essentially consistent with zero for 
L5 > 16. We expect the shifts in the nucleon mass induced by these small shifts in the pion mass to be negligible, 
and indeed. Figs. ^ and ([7]) show that the differences between the nucleon mass at a general value of L5 and at 
L5 = 16 are consistent with zero for L5 > 16. Hence, we choose L5 = 16 to be a good compromise between accuracy 
and performance. 
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FIG. 3: Residual quark mass as a function of Ls for the two samples (heavy and light) of 25 configurations each. 
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FIG. 4: Dependence of the pion mass on the extent 
of the fifth dimension L5 for heavy quarks. 



FIG. 5: Dependence of the pion mass on the extent 
of the fifth dimension Ls for light quarks. 



B. Tuning the quark mass 

We define the light quark masses in our hybrid theory by matching the pion mass in two calculations: (i) using two 
plus one flavors of dynamical Asqtad sea fermions and Asqtad valence fermions|29j and (ii) using the pion mass in 
our hybrid calculation with Asqtad dynamical sea fermions and valence domain-wall fermions with L5 = 16. Because 
of the four tastes and correspondingly sixteen light pseudoscalar mesons in the staggered theory, it is necessary 
to choose between matching the lightest pseudoscalar mass, corresponding to the Goldstone pion of the theory, or 
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FIG. 6: Dependence of the nucleon mass on the FIG. 7: Dependence of the nucleon mass on the 

extent of the fifth dimension L5 for heavy quarks. extent of the fifth dimension L5 for light quarks. 
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dataset 




# 




(am)g 






/ \Asqtad 

(a"^)]v 


(am)iv 


m°^P [MeV] 


1 


20"* X 32 


425 


0.050/0.050 


0.0810 


0.4836(2) 


0.4773(9) 


1.057(5) 


0.986(5) 


758.9(1.4) 


2 




350 


0.040/0.050 


0.0478 


0.4340(3) 


0.4293(10) 


1.003(3) 


0.938(8) 


682.6(1.6) 


3 




564 


0.030/0.050 


0.0644 


0.3774(2) 


0.3747(10) 


0.930(3) 


0.869(6) 


595.8(1.6) 


4 




486 


0.020/0.050 


0.0313 


0.3109(2) 


0.3121(11) 


0.854(3) 


0.814(7) 


496.2(1.7) 


5 




655 


0.010/0.050 


0.0138 


0.2242(2) 


0.2243(10) 


0.779(6) 


0.730(12) 


356.6(1.6) 


6 


28^ X 32 


270 


0.010/0.050 


0.0138 




0.2220(9) 




0.766(15) 


352.3(1.4) 



TABLE I: The lattice volume Q and number of configurations used for the DWF calculations and a comparison of the quark, 
pion, and nucleon masses in the DWF and Asqtad calculations as described in the text. 

some appropriately defined average. In this work, we have chosen to match the Goldstone pion, and the results of 
tuning the domain wall quark mass such that the domain wall pion mass agrees within one percent with the Asqtad 
Goldstone pion mass are shown in Table HI The substantial difference between the bare quark masses for Asqtad and 
DWF valence quarks reflects the significant difference in renormalization for the two actions. An observable physical 
difference is the fact that once the DWF quark masses have been adjusted to fit the Asqtad Goldstone pion masses, 
the DWF nucleon masses are approximately 6% lower than the corresponding Asqtad nucleon masses. We attribute 
this to the range of pseudoscalar masses in the staggered theory and note that had we used a heavier quark mass 
so that the DWF pion fit some appropriately weighted average of the staggered pion masses, then the DWF nucleon 
would have been heavier. 



C. Operator renormalization 

The quark bilinear operators in Eq. [5] are renormalized using a combination of one-loop perturbation theory and 
non-perturbative renormalization of the axial vector current. 

Our lattice calculations using lattice regularization with cutoff 1/a are related to physical observables at scale /i^ 
in the MS renormalization scheme in 1-loop perturbation theory by 

= ■ 0^^^^{a% (17) 

where the anomalous dimensions and the finite constants Bij have been calculated for domain wall fermions with 

HYP smearing in Refs. [3^, |4^. Because the renormalization factors for operators with and without 75 are identical 

at quark mass zero, we use mass independent renormalization with all renormalization constants defined at quark 

mass zero. The renormalization factors ZP, for domain wall mass M = 1.7 used in this work are tabulated in Table 
*j 

im using the results for the one-loop coupling constant g^/(127r^) — 1/53.64 from Refs. [33, |40|. By virtue of the 



operator 


H{4) 


^0,pcrt 


Qbshlt^D^Dpyq 
Qbshit^D^Dpyq 




0.962 
0.968 
0.980 
0.982 



TABLE IL Perturbative 1-loop lattice renormalization constants for the MS scheme at a scale /i^ = 1/a^. 

suppression of loop integrals by HYP smearing, the ratio of the one-loop perturbative renormalization factor for a 
general bilinear operator to the renormalization factor for the axial current is within a few percent of unity, suggesting 
adequate convergence for this ratio at one-loop level. Since one element in the calculation common to all operators 
arising from the wave function renormalization in the fifth dimension is not small, it is desirable to determine this 
one common factor non-perturbatively. This is accomplished using the fact that the renormalization factor, Za, for 
the four dimensional axial current operator — q^g^^a rnay be calculated using the five dimensional conserved 
axial current for domain wall fermions by the relation [36l| {A^{t)A^{Q)) = Z A{A^{t)A^{{))) . Hence the complete 
renormalization factor is written as the exact axial current renormalization factor times the ratio of the perturbative 
renormalization factor for the desired operator divided by the perturbative renormalization factor for the axial current. 
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FIG. 8: Unpolarized (vector) generalized n = 2 form factors for the flavor combinations u — d (left) and u + d (right). 
Disconnected contributions are not included. 
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In the continuum, because of Lorentz invariance, the totally symmetric operator p^q cannot mix with 

the mixed symmetry operator g[75]7[^D{^]Z?p}.(;, where the square brackets denote antisymmetrization. In contrast, 

on the lattice, both operators appear in the same representation, r^^ ^ so that they can and do mix. However, the 



mixmff coefficient [31140], Zf^ = 2.88 x lO'^, is very small, so that we have ignored the contribution of the mixed 
symmetry operator in this present work. 

All results below have been transformed to a scale of //^ = 4 GeV^ 
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FIG. 9: Polarized (axial vector) generalized n = 2 form factors for the flavor combinations u — d (left) and u + d (right). 
Disconnected contributions are not included. 



IV. NUMERICAL RESULTS FOR THE GENERALIZED FORM FACTORS 

Since two point functions taken at the sink r — Tgnk, C'^^'^i'Tsnk, P') and C^P*(rsnk, -P) in the ratio Eq. (fT5|) decay 
exponentially for the full Euclidean distance Tgnk — Tsrc , they are particularly subject to statistical noise. In the 
worst case, they may become negative, which we observe for three values of the momentum transfer for the dataset 
m = 0.01, 20^^. The corresponding datapoints are excluded from our analysis. Our numerical results for the complete 
set of unpolarized and polarized n = 1,2,3 isovector and isosinglet GFFs as functions of the momentum transfer 
squared are provided in appendix [XI 
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FIG. 10: Flattening of the slope of the A„o GFFs with increasing n for flavor combinations u — d (left) and u + d (right). The 
solid curves and error bands correspond to dipole fits described in the text. Disconnected contributions are not included. 

In Figs. ([5]) and ^ we show results for the vector generahzed form factors j42o, S20, C20 and axial vector GFFs 
^20, B2Q as functions of the momentum transfer squared t. We observe that the absolute values in the isovector and 
isosinglet channels are in qualitative agreement with the predictions from large Nc counting rules, see e.g. [4l[, for 
the unpolarizcd GFFs 
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In the polarized case, the inequalities from the counting rules are not satisfied nearly as strongly. Whereas the counting 
rules predict: 
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FIG. 12: Two dimensional rms radii of the axial 
vector GFFs versus for the flavor combination 
u-d. The results for = 354 MeV, = 20^ are 
displayed in gray and slightly shifted to the right 
for clarity. 



our results for O'^ly slightly larger than ^20^'' > ^'^d although the errors are large, appears to be 

comparable to i?2o'^'' rather than dominating it. Finally, our results disagree with the predicted hierarchy between 
different types of GFFs: 

iSaV'l - » \A-^+^\ ^ Nl , (21) 

since the lattice results (at non-zero t) clearly give ^20^'' > B2q'^. It would be valuable to understand why these 
counting rules are only partially satisfied. 

For future reference, it is important to note that the GFF C20 , which gives rise to the ^-dependence of the n = 2 
moment of the GPDs H{x,£,,t) and E{x,£,,t), is compatible with zero for u — d, over the full range of momentum 
transfer squared t w —0.12 ... — 1.2 GeV^. Similarly, the isosing let GFF 5^0+'*, which is one of the terms in the 
contribution of the total angular momentum to the nucleon spin, is compatible with zero within errors. We will study 
both these GFFs in detail in section V. 

We now consider the behavior of the slopes of the GFFs Ano and their relation to the transverse size of the nucleon. 
Since J^^ dxx'^~^H{x, ^ = 0, t) = Ano{t), it is evident that the GFFs for increasing n correspond to increasing average 

momentum fractions (x). As the average momentum fraction gets larger, or equivalently as n ^ 00, we expect the 
t-slope of the GFFs Ano to flatten. This may be understood in terms of the light cone Fock representation [i^, by 
the fact that the final state nucleon wavefunction for a struck quark with momentum fraction x and initial transverse 
momentum fc"' depends on the transverse momentum k± = A:™ — (1 — x)A±. Hence, a large transverse momentum 
transfer t = — can be better absorbed without causing breakup of the bound state by quarks with large momentum 
fraction x. Additional insight is obtained by considering the impact parameter dependent GPD, H(x,b\), which 
has a probability interpretation and is the Fourier transform Q with respect to transverse momentum transfer of 
i/(x,e = 0,t = -Ai): 

H{x, hi) = j ^e-^'^-^-Hix, ^ = 0, i = -Ai) , (22) 

where A_l is the transverse momentum transfer. The impact parameter b± corresponds to the distance of the active 
quark from the center of momentum of the nucleon. As a; ^ 1, a single quark will carry all the longitudinal momentum 
of the nucleon and therefore represent its center of momentum, so that the impact parameter distribution in this limit 
is strongly peaked around the origin, H{x — > 1,65_) '^^(^D- The corresponding flattening of the GFFs in the 
momentum transfer t is clearly visible in Fig. (jlOp . where we compare the slopes of the GFFs ^(n=i,2,3)o which have 
been normalized to unity at t = 0. 

Dipole fits to the GFFs in Fig. (fTU)) . denoted by the solid lines and statistical error bands, enable us to determine 
the slopes of the form factors and thus express the two- and three-dimensional rms radii (r^)^^ and (r^)^/^ in terms 
of the dipole masses niD 
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FIG. 13: Ratio of generalized form factors y43o(t)/Aio(t) for the flavor combinations u 
with the parametrization in Ref. [lj|. Disconnected contributions are not included. 



d (left) and u + d (right) compared 



Since the range of values for the momentum transfer t is much smaller for the large volume {L^ — 28^) dataset, we 
have restricted the dipole fits for all datasets to the overlapping region oi t ~ . . . — 0.8 GeV^. Our results for the 
2d rms radii versus the pion mass squared are presented in Figs, pip and p2p . These results confirm the dramatic 
dependence of the transverse rms radius on the moment and thus the average momentum fraction as first observed [l3| 
for pion masses 750 MeV and higher, and show that this dependence increases as the pion mass decreases. Indeed, 
considering the ratio of the n = 3 moment to the n — 1 moment, which both correspond to the same sum or difference 
of quarks and antiquarks, we observe that for vector GFFs this ratio decreases from approximately 0.58 to 0.22 as the 
pion mass decreases from 750 MeV to 350 MeV, and for axial vector GFFs, it decreases from roughly 0.71 to 0.43. 

In Figs. and ^1)) we present a first comparison of our results for ratios of generalized form factors A3o{t)/Aio{t) 
and A3Q{t) / Aio{t) to the parametrization by Diehl et al.[Tl| as function of the momentum transfer squared t. As the 
pion mass decreases, the slope of our results approaches that of the phenomenological parametrization. Our results 
clearly indicate that a factorized ansatz for the GPDs in x and t, which would lead to constant ratios in Figs. (fTT 
and (jl4p breaks down already for small values of the momentum transfer squared |t| ^ 1 GeV^. 
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FIG. 14: Ratio of polarized generalized form factors yl3o(t)/Aio(f) for the flavor combinations u — d (left) and u -\- d (right) 
compared with the parametrization in Ref. J4|. Disconnected contributions are not included. 



The GFFs A^^it = 0) = {x)q and i?|o(^ ~ 0) enable us to compute the total quark angular momentum contribution 
to the nucleon spin [ll, J'' = 1/2{A\q{{)) -|- i?2o(0)). Figures (fTS)) and HH) show results for the quark spin A\Q{t = 
0)/2 = AYfl /2 and the orbital angular momentum L'^ = J'^ — contributions to the nucleon spin S — 1/2 

versus the pion mass squared. Preliminary chiral extrapolations of AT,'' based on self-consistently improved one-loop 
ChPT [4i,l4i,|4i,[i3,|4i] for AE"+'' and ChPT including the A resonance [H, [53| for gA = AE""^ and are shown 
as shaded bands. The values for i?|o(^ = 0) have been obtained from a linear extrapolation of B^^'^'it) and a dipole 
extrapolation for -620'^ (i) in t. The resulting uncertainty in i3|o(* ~ 0), which contributes to the uncertainty in L^, 
depends on the details of the corresponding fit, such as the functional form and range of i, and is therefore partially 
systematic. To allow the reader to assess the absolute statistical errors, we represent the errors for Lq coming from 
the extrapolation in t by error bands around the m^-axis in Figs. p3|) and Experimental results for the quark 

spin fractions AE"+'' and A(w, d) — AE"''' are represented by open stars for the prediction given in the HERMES 
publication from 1999 ^5lj and filled stars for the 2007 HERMES results [s^l- The significant difference between the 
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new HERMES results, which are consistent with recent COMPASS results[53, and the values given in 51| is probably 
to a large extent due to the simple Regge-parametrization which has been used in (sij to compute the contribution 
to AE coming from the low x- region. It is gratifying that the new values are much closer to our lattice results. 

These results reveal two remarkable features of the quark contributions to the nucleon spin. The first is that the 
magnitude of the orbital angular momentum contributions of the up and down quarks, and L"^, are separately 
quite substantial, starting at 0.15 at m^r = 750 MeV and increasing to nearly 0.20 at 350 MeV, and yet they cancel 
nearly completely at all pion masses. The second is the close cancellation between the orbital and spin contributions 
of the d quarks, and for all pion masses. It would be valuable to understand the physical origin of both 

features. 
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FIG. 15: Total quark spin and orbital angular momentum 
contributions to the spin of the nucleon. The filled and 
open stars represent values given in HERMES 2007 [s^ and 
1999 [sH ] respectively and open symbols represent earlier 
LHPC/SESAM calculations. The error bands are explained 
in the text. Disconnected contributions are not included. 



FIG. 16; Quark spin and orbital angular momentum contri- 
butions to the spin of the nucleon for up and down quarks. 
Squares and triangles denote AE" and AS'* respectively, and 
diamonds and circles denote L" and L"^ respectively. The 
filled and open stars represent values given in HERMES 2007 
[S^ l and 1999 ^Sjl respectively and open symbols represent 
earlier LHPC/SESAM calculations. The error bands are ex- 
plained in the text. Disconnected contributions are not in- 
cluded. 



V. CHIRAL EXTRAPOLATIONS 



Our ultimate goal is to use the combination of full QCD lattice calculations in the chiral regime and chiral pertur- 
bation theory to extrapolate to the physical pion mass, to extrapolate to infinite volume, to extrapolate in momentum 
transfer, and to correct for lattice artifacts, with all the relevant low energy constants being determined solely from 
lattice data. Significant progress has been made in many aspects of chiral perturbation theory (ChPT) relevant to 
the nucleon observables addressed in this work [i^, HO, 13, [H, HI, 113, [H, [H, 6C , 61, 62, "Mj [H, Hg]. Although 
important developments have been made in correcting for our hybrid action [61 



62, 63, 
62, M, M, 




and finite volume 



[iO, [S^i results for the relevant GFFs are not yet available. In this work we will focus on ChPT treatment of the pion 
mass and momentum dependence. 

The basic problem is that currently, there is not yet unambiguous evidence supporting a particular counting scheme 
and its convergence criteria, leading to a range of alternative re-summations, and there is similar ambiguity concerning 
the choice of degrees of freedom, such as when and if it is essential to include the A resonance. When complete results 
for the observables of interest are available, it will be interesting to compare four approaches: heavy baryon ChPT 
(HBChPTlE" 



Si EE 




, covariant ChPT in the baryon sector (BChPT ')[65i , self-consistently improved one- loop ChPT 
and ChPT with finite-range regulators|57l Isol |66| |. Although self-consistent improvement by 
utilizing values of parameters like /jr and qa calculated on the lattice at the relevant pion mass and finite-range 
regulators appear to improve the behavior of ChPT at larger values of the pion mass, based on the results available 
in the literature, we will focus on the two formulations HBChPT and BChPT. 

Heavy baryon ChPT, which wc will subsequently always refer to as HBChPT, assumes that to^ and the magnitude 



of the spatial three- momentum, p, are much smaller than the nucleon mass in the chiral limit. 



'■N 



890 MeV, 



and the chiral scale A^ = '^t^J-k — 1-17 GeV, and simultaneously expands in powers of the four quantities e = 
{yf~' X^' contrast, covariant baryon ChPT, which, slightly changing the notation of Ref. [6^, we will 

subsequently always refer to as CBChPT, does not treat and as comparable scales, but rather keeps all powers 
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f)i [GeV] 


m'l^ [GeV] 






ci[GeV-'] 


c2[GeV-'] 


C3[GeV-^] 




A = m'A - m% [GeV] 


0.092 


0.89 


1.2 


0.5 


-0.9 


3.2 


-3.4 




0.3 



TABLE III: Low energy constants used in the chiral extrapolations. 



(^3-)" generated by the couplings included in the ChPT Lagrangian. Thus, it is a resummation that includes terms 

that would contribute in higher order to HBChPT and may be thought of as recoil corrections. The HBChPT results of 
refs. ^62, 63, 64] and the CBChPT results of ref. [6^ have the desirable property that they use the same regularization 
scheme, so that truncation of the higher order terms in CBChPT yields the corresponding result in HBChPT, a 
feature that we will utilize below. One of our primary objectives will be to assess the regimes of applicability of these 
two alternative formulations. For notational convenience, we will refer to the generic momentum in both theories 
as p, and order both theories in powers of p". We would like to note, however, that the counting scheme of the 
HBChPT-approaches in l6p . [63l. [64] differs from the one used in the CBChPT-approach of [65| . 

The HBChPT [i^, m, and CBChPT [6^ results for GFFs, including the dependence on the momentum transfer 
squared, t, enable us to investigate for the first time possible non-trivial correlations in the TOtt- and i-dependence of 
GFFs. It is interesting to note that to order 0{p^) in HBChPT, unpolarized and polarized GFFs are independent 
of each other and depend on separate chiral limit values and counter-terms. In contrast, in CBChPT, the pion-mass 
dependence of the isovector momentum fraction of quarks {x)u-d is simultaneously controlled by both the chiral 
limit values {x)^_^ and (a;)^„_^^, as is also the case for {x)Au-Ad- To 0{p'^), however, CBChPT does not include 
insertions from pion operators, and it turns out that the i-dependence for the isosinglet (and isovector) case to this 
order is therefore essentially linear and decouples from the pion mass dependence. Once CBChPT calculations have 
been pushed to higher orders, it will be interesting to study the combined non-analytic (t, m^)-dependence of the 
full set of polarized and unpolarized GFFs based on this approach. For the time being, we will investigate possible 
non-trivial correlations in TOtt and t in the framework of covariant CBChPT by including partial C'(p'^)-corrections as 
discussed below. 

The low energy constants used for the chiral extrapolations are summarized in Table [Till Ultimately, all these values 
will be determined simultaneously by a global fit to a full set of lattice calculations of all the relevant observables using 
the same lattice action, but at present they are chosen as follows. We will use /° — 0.092 GeV, m% — 0.89 GeV and 
g'^ = 1.2 as the chiral limit values of the pion decay constant, nucleon mass and the axial vector coupling constant, 
respectively, and for notational simplicity, will subsequently omit the upper index for these quantities. We note that 
these values are compatible within statistical errors with hybrid lattice calculations of rriN, and qa- In addition, 
we need the low energy constants ci, C2 and C3 for the chiral extrapolation of Cjo^'^ in the framework of HBChPT 
in section IVl Here, cs = 3.2 GeV'^ and C3 = -3.4 GeV-^ are taken from Refs. [fii, [tiI and a = -0.90 GeV-^ 
has been obtained from a CBChPT fit to our nucleon mass lattice data, which provides us with a parametrization 
for the pion mass dependence of the nucleon mass in our simulation needed for some of the chiral extrapolations 
below. Depending on the order of ChPT, diagrams with insertions of pion operators contribute for the isosinglet 
GFFs, which introduces the momentum fraction of quarks in the pion in the chiral limit, {x)'^'^^, as an additional low 
energy constant. Ultimately, we will calculate this quantity from chiral fits to hybrid lattice results for the pion, but 
for now we use (a;)^^j = 0.5 [t^. ItsI . . which is obtained from lattice calculations that are in reasonable agreement 
with phenomenology ' 75|, Wdi. l77l] . Finally, for the chiral extrapolation of the total quark angular momentum in the 
framework of HBChPT including the A-resonance, we use the nucleon-A mass splitting A — 0.3 GeV, which is 
consistent with lattice calculations and the large-A^c relation = 3/{2-^/'^)g^ for the pion-nucleon-A couphng 

dirNA- This latter result will soon be superseded by extrapolation of lattice calculations of the iV — A transition form 
factor HlliiOl. 

Our chiral extrapolations are organized as follows, and summarized in Table IIVI We begin in section IV Al with a 
comparison of CBChPT and HBChPT extrapolations of the isovector GFF A2Q'^{t) and show that whereas CBChPT 
yields a satisfactory fit over the range of pion masses used in the lattice calculations, HBChPT only produces fits to the 
lowest few points. Hence, in sections fV BI through lV F[ we study the pion mass and i-dependence of the isovector GFFs 
BI^q'^ and C^^'^ and of the isoscalar GFFs ^20^''' ^20"^ ^.nd C^^'^ . This is followed in scction lVTJl bv a discussion of our 
results for the angular momentum of quarks, based on the CBChPT extrapolations for ^20 = 0) and i?2o(* = 0). In 
the counting scheme of [gl, HI, HI] insertions of pion operators occur at OijP') in HBChPT, leading to a non-analytic 
combined dependence on m.^- and t for the GFFs B!^+'^{t) and G^^'^{t) , which we study in sections |VHl [Vl andlVJ] 
Finally, in section IV Ki we study the pion mass dependence of the total quark angular momentum J, in HBChPT, 
both including and excluding explicit A degrees of freedom and compare with the corresponding CBChPT results. 
The chiral extrapolations of B'!^+'^{t) and G^^'^it) in sections EB EH andlVJl are the first parametrizations of their 
combined (m^, t)-dependence, providing valuable insights into non-trivial correlations of m,r and t. Although we have 
a considerable amount of lattice data available for reasonably small values of \t\ < 0.25GeV^ and < SOOMeV, we 
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section 


GFF 


HBChPT 


CBChPT 


expected dependence on m.^,t 


A 


AU — d 

^20 






non-analytic in m^, « linear in t 


B 


j-yU — d 

^20 




0(p^) + corr. of ©(p^) 


non-analytic in m^, ~ linear in t 


C 


/^ru — d 
'^20 




0{p^) + corr. of ©(p^) 


non-analytic in m^, ~ linear in t 


D 


^20 




0(p^) + corr. of ©(p^) 


non-analytic in m-^ and t 


I, E 


r)u+d 
^20 




O(p^) + C'(p^)-CTs 


non-analytic in m-n and t 


F 


f^u+d 
•-^20 




O(p^) + corr. of ©(p^'"*) 


non-analytic in m-^ and t 


G 


ju+d ^ 1/2(^ + 5)^0+" 




O(p') + corr. of ©(p^) 




H 


E^+'' = {A + t/(4mMrB)^+^ 


0{p') 




linear in and t 


H 








non-analytic in and t 


J 


/-iu-\-d 
'-'20 


O(p') 




non-analytic in ra-n and t 


K 


ju+d ^ 1/2(^ + 5)^0+" 


O(p') 






K 


ju+d ^ i/2(A + B)^o^" 


0{p^) with A 







TABLE IV: Overview of different approaches to the (m^, t)-dependence of GFFs in ChPT studied in sections IV Al V Kl 



use an extended set of results for \t\ < 0.48GeV^, TOtt < 700MeV in most of the fits to improve the statistics. 

A. CBChPT extrapolation of A^o""(t) 
The 0(j)^) CBChPT result [H for the isovector GFF A'i^^'^it) is 

where /^"''(w^), "v) and j^'^im.,^) contain the non-analytic dependence on the pion mass and momentum 

transfer squared and ^20" = ~ 0, ttItt = 0). Because of the small prefactor, the term cx hA{t,niT^) is 

of O(10~^) for rn^ < 700 MeV, |t| < 1 GeV'^ and therefore numerically neghgible. Thus, there are essentially no 
correlations of t and TOtt present, and the dependence on t is only due to the counter term (^20'"'* We use the value 
^20"""^ = 0.17 obtained from a chiral fit to our lattice results for Ai^g'^it = 0) = {x)Au-Ad Since the low energy 
constant A^q""'^ is a common parameter in the CBChPT-formulae for the GFFs ^20''^' -^20'^ ^^'^ ^2^'^^ performed 
a simultaneous fit based on Eq. (EH) (for A^^'^), Eq. ^ (for B^^'^) and Eq. ^ (for C^q'^) with a total of 9 (1 
common and 8 separate) fit parameters to over 120 lattice datapoints. The details of the CBChPT-extrapolations 
and the results for the GFFs -62(7^ ^'^d C20 will be discussed below in sections IV Bl and IV Ci respectively. We 
find A20" = 0.133(9) and {x)u-d = A'^^'^it = 0, m7r,phys) = 0.157(10) at the physical point. This is in very good 
agreement with phenomenological results from CTEQ and MRST ^1] PDF-parametrizations, (a;)jfj^^'^^°°^ = 0.157(5) 

and {x)1T^^^ = 0.155(5). The results of the fit are shown in Figs. ^ and The dependence of Al^'^{t) on 

the momentum transfer squared is presented in Figs. (fTO|) and (I^Hl) . where we again obtain a good description of the 
lattice data. 

To study the difference between HBChPT and CBChPT, we took the heavy baryon limit of CBChPT while 
keeping the same values of the fit parameters, and obtained the dotted line in Fig. (I19p . This curve only overlaps the 
CBChPT curve for < ?7i.7r,phys and drops off sharply for > m7r,phys, indicating the quantitative importance 
of the truncated terms when using the coefScients from the CBChPT fit. In addition, it is important to ask the 
separate question of how well the lattice data can be fit with the HBChPT expression when the coefficients are 
determined directly by a best fit to the data. The dashed curve in Fig. p7|) shows the result of fitting our lattice data 
for \t\ < 0.3GeV^ and 171^^ < 0.5GeV, and indicates that HBChPT describes the behavior of our lattice data over a 
significantly smaller range of pion masses than CBChPT. 

Because limitations in computational resources presently require us to include lattice data extending to such large 
pion masses, it would clearly be desirable to carry out a chiral perturbation theory analysis consistently including 
all terms of C'(p^). In the absence of the requisite full ChPT analysis, we have studied uncertainties in the chiral 
extrapolations by repeating the fit for different maximal values of the included pion masses. Figure (j2ip shows a 
comparison of the chiral extrapolations of A^q''', based on fits to the lattice data in the regions TOtt < 500, 600, 685 
and 760 MeV. We find that the extrapolations to the chiral limit fully agree within statistical errors in all four cases. 
Note that the experimental point shown in Fig. ((2T|) was not included in the fits, but each of the four analyses is 
consistent with it. This insensitivity to the upper mass cutoff shows that the strong bending towards the physical 
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FIG. 



17: Lattice results for A^q"^ at i = GeV^ 



versus . The error band is the result of a global 
simultaneous chiral fit using Eqs. (|24|l . (I25|l and 
(|26|) . The phenomenological result from CTEQ6 is 
indicated by the star. The heavy-baryon-limit of 
the CBChPT fit is shown by the dotted line, and a 
HBChPT fit to the lattice data for \t\ < O.SGeV^ 
and m-n < 0.5GeV is shown by the dashed line. 



FIG. 18: Lattice results for A^q at t « -0.24 
GeV^ versus together with the result of a 
global simultaneous chiral fit using Eqs. ((24]), (I25p 
and flE} . 
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FIG. 19: Lattice results for A2o'^{t) at « 0.35 
GeV together with the result of a global simulta- 
neous chiral fit using Eqs. ((25)l and ((26)l . 



FIG. 20: Lattice results for A2o''(t) at « 0.496 
GeV together with the result of a global simulta- 
neous chiral fit using Eqs. ((24)) . ((25j and ((26j. 



point is not driven by the large pion mass region, where 0{p'^) corrections would be largest. Furthermore, all four 
chiral fits are, within statistical errors, in agreement with the lattice data points at large pion masses. This indicates 
that the present statistical error envelope is comparable to any systematic effects due to higher order corrections. 

Another prescription to estimate 0{p^) corrections that has been advocated in the literature Ref. [1^ is simply 
adding a single term and, assuming both "naturalness" of the coefficient and the lack of other functional forms, 
seeing what error band arises from varying the coefficient from -1 to -1-1. Thus, to explore this possibility, following 
Ref. [sHi, we have added to the result in Eq. ((Mj) the term (5^'''"~''m^/(A^TO^), where A^^ « 1.2 GeV is the chiral 
symmetry breaking scale, and varied the constant (J^'*'" in the range — 1, . . . , +1. The results of fits to the lattice data 
for A2q'^ including this additional term are shown in Fig. ((Hj), where the error band corresponds to (5^'''" = 0, the 

dashed line corresponds to (5^ = +1, and the dotted line corresponds to = —1. From the figure, we note 

that this term alone with coefficients +1 and -1 is clearly inconsistent with the behavior of the data. Theoretically, 
this is not unreasonable, since the foundation of the IR regularization scheme is a resummation of classes of terms, 
and here a single cubic term has been arbitrarily singled out. Analogous fits with similar qualitatively inconsistent 
behavior were also obtained for A2^'^, treated in a later section. Hence, although this prescription may provide useful 
estimates in other contexts, we do not believe it is useful in this work, and hence we do not include it in subsequent 
fits. 
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FIG. 21: Lattice results for A^^'^ at i = GeV^ 
versus together with chiral fits based on 

Eq. (|24|l . The four different error bands represent 
chiral fits to lattice results including pion masses 
in the regions < 500, 600, 685 and 760 MeV. 
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FIG. 23: Lattice resuhs for B^^'' at t « -0.24 
GeV'^ versus together with the result of a global 
simultaneous chiral fit using Eqs. (I24|l . (|25|l and 
(ESI). 
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FIG. 22: Chiral fits including the terms in Eq. pi|l 
plus an additional term as described in the text. 
The narrow band is the original band omitting this 
term, the dashed line corresponds to 5^^^'^~'' — 1, 
and the dotted line to 5^-''""'' = —1. 



0.1 









m,=.0.35GeV 






I'l 





0.2 0.4 0.6 0.8 1 1.2 

-t [GeV^] 



FIG. 24: Lattice resuhs for B^g''{t) at k: 350 
MeV versus (— t) together with the result of a global 
simultaneous chiral fit using Eqs. (|24p . (|25p and 
(EH. 



B. CBChPT extrapolation of B^^'^{t) 

The CBChPT calculation ^ for the isovector B20 GFF gives 

B^^\t, m,) = + "^-) + "^^("^-) |^^"-'^^ + S^-'^-^rrA , (25) 

TTLN ruN I J 

where mAr(m7r) is the pion mass dependent nucleon mass, B%^~'^ = B^^'^{t = 0,TO^ = 0), and where we have 
included estimates of O(p^)-corrections in form of {5^^~'^ t) and (5^'"'""'* m^). The low energy constants i?2o"~'', 
^^n-d grnT,M-d treated as free parameters and may be obtained from a fit to the lattice data. The non-analytic 
dependence on m^r and t is given by m,r)j but it turns out that this function is approximately independent 

of t, /i^~''(t,m^) K, /i^"''(m7r) for < 700 MeV, \t\ < 1 GeV^. The t-dependence is therefore in practice Hnear 
due to the C'(p'^)-correction term. For the WTr-dependent nucleon mass we use 0{p'^) CBChPT [6l,[84| fitted to our 
lattice results for mjv- The chiral extrapolation of ^20''^ based on a global simultaneous fit as discussed in the 
previous section with A^q""'' as common fit parameter in Eqs. ([M]) . ([^5]) and (|26p . We obtain i?2o"~'^ = 0.263(62) 
and B^^'^it = 0, m^^phys) = 0.273(63) at the physical pion mass. Results of the fit are shown in Figs. (|^ and ([M)) . 
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C. CBChPT extrapolation of C^^'^it) 



The pion mass dependence of the isovector GFF C20 in CBChPT to 0{p^) is very similar to that of the isovector 
B20 above and given by [6^ 

C20 (t,m^) = 0*20 h(j [t,m^)-\ i6(^ t + rn^\ , (26) 

m,M rriM I J 

where 6*20""'^ = C2Q'^{t — 0,to^ — 0). As in the case of i?20 ((5^"^''^) and ((5™'""^'' m^) represent C'(p'^)-correction 




FIG. 25: Lattice results for C^o at t fa -0.24 
GeV^ versus together with the result of a global 
simultaneous chiral fit using Eqs. (|24|l . (|25|l and 



FIG. 26; Lattice results for C2o''{t) at f« 350 
MeV versus {—t) together with the result a global 
simultaneous chiral fit using Eqs. (P^ . and 
If26ll. 



terms, and it turns out that h^'^ (t, m^) is practically independent of t. From a global simultaneous chiral fit based on 
Eqs. and ([^5]) with common fit parameter ^20'"'* discussed in section lV A[ we obtain C20" = —0.017(39) 

and C2o~''(i = 0,TO7r,phys) = —0.017(41) at the physical point. The results are presented in Figs. (^5]) and ^E^. We 
note that C^o^^it) is roughly one order of magnitude smaller than A^q '^(t) and '^{t) and fully compatible with zero 
within errors, C20~'^(*) ~ 0. This implies a rather mild dependence of the n = 2 moment of the GPDs iJ"~'^(x, ^, i) 
and E"^'^{x,S^,t) on the longitudinal momentum transfer ^ (at least for small ^), so that H'^^^{S^,t) w ^2(7'^('') 



D. CBChPT extrapolation of A^+'^{t) 

The (total) isosinglet momentum fraction of quarks, Al^g'^it = 0) = {x)u^ii is not only an important hadron structure 
observable on its own but is in addition an essential ingredient for the computation of the total angular momentum 
contribution of quarks to the nucleon spin, J"+'' = l/2{A2a'^i^) + -S20^'^(0))- The combined (i, m7r)-dependence in 
CBChPT is given by [H: 

Al^+'^it, m.) = A%-+' - ^^)) + ^20"'"+'™^ + 4o'+'^ + ^Al^tit, rn^) + 0{p') , (27) 

where ^20"^'' = ^'io'^i^ ~ 0,m^ = 0), /^^'^(m-ir) and hA{t,mTr) contain the non-analytic dependence on the pion 
mass and momentum transfer squared, and the constants A^q"'"'^'^ and y4*'"+'^ may be obtained from a fit to the 
lattice data. In this counting scheme, contributions from operator insertions in the pion line proportional to the 
momentum fraction of quarks in the pion in the chiral limit, (x)^'^^, are of order 0(j)^). However, in order to see 
if such contributions could be relevant for the pion masses and values of the momentum transfer squared accessible 
in our calculation, we include the estimate of the C'(p'^)-contribution AAI^q"^ provided in [6^ in the fit to the lattice 
data points. 

Similar to the isovector case discussed in the previous sections, the low energy constant A'^'q'^'^ is a common 
parameter in the chiral extrapolation formulae for the isosinglet GFFs ^20''' i ^20 ^^'^ ^20'^- Using (a;)^^^^ = 0.5 
from Table UTTl as an input parameter, we performed a simultaneous fit to over 120 lattice data points for these three 
GFFs, based on Eqs. ([?7|) . ([^5]) and ([^5]) . with 1 common and 8 separate low energy constants as fit parameters. For 
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the details of the CBChPT results for B^^ 



A^o""*"* = 0.524(25) and 



and 



we refer to the sections IV El and IV Fl below. The chiral fit gives 

x)v.+d — A^n'^{t — 0, m7r,phys) = 0.520(24) at the physical point. Again, this is in very good 

0.538(22) 



agreement with phenomenological results from CTEQ and MRST [81[ parametrizations, (2;)^^" 



MRST2001 



and (x 



CTEQ6 



: 0.537(22). A variation of the input parameter (x)^!^^ by ±10% only leads to a small change in 

' 0(1%), which is significantly smaller than the statistical error of w 5%. The results of the fit are 
27l) and ((28)1 . We would like to note that the slight upwards bending in Fig. ([27l) at low ttIt^, and 



shown in Figs 

therefore the good agreement with the phenomenological value, is due to the C'(p'^)-contribution AA2q'^. It has to be 
seen if this somewhat unusual curvature persists once the full 0{p^) contribution is available and fitted to the lattice 
results. The inclusion of contributions from disconnected diagrams could also require a different extrapolation in 771,^. 
The dependence of A2Q'^{t) on t at fixed values of m^r is presented in Figs. and pop . 

As in the case of A^^'^, we also consider the heavy baryon limit of the CBChPT fit, giving the result A'l^^'^it = 



0,m^) = A%^+'' + A 



4 .u-\-d 
^20 



represented by the dotted line in Fig. l|77|) . which agrees with the CBChPT result only 



over a very limited range at low pion masses. Notably, while the lattice results for are rising for larger pion 



masses, the heavy-baryon-limit curve has the opposite slope with negative A^"'"^''. 



However, a direct HBChPT fit 



with free coefficients (see also section IV H|) shown by the dashed curve leads to a positive 
description of the lattice datapoints. 



and a reasonable 
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FIG. 27: Lattice results for y4^+'* at t = GeV^ 
versus m\. The error band is the result of a global 
simultaneous chiral fit using Eqs. (|27|l . (|28|l and 
(|29p . The phenomenological value from CTEQ6 is 
denoted by a star. The heavy-baryon-limit of the 
CBChPT fit is shown by the dotted line, and a 
HBChPT fit to the lattice data for \t\ < O.SGeV^ 
and < 0.5GeV is shown by the dashed line. 



FIG. 



"^{t) at t « -0.24 
the result of a 
global simultaneous chiral fit using Eqs. (|27p . (I28|l 
and (p9)) . 



28: Lattice results for 
GeV'^ versus together with 
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FIG. 29: Lattice results for Al^'^{t) at 



0.35 



GeV together with the result of a global simulta- 
neous chiral fit using Eqs. ([27l, ((28l) and ((29)l . 
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FIG. 30: Lattice resuhs for A^^'^{t) at 



; 0.496 



GeV together with the result of a global simulta- 
neous chiral fit using Eqs. ((27l), JSH and (|29|l . 



22 



E. CBChPT extrapolation of B^+''{t) 



The dependence on the pion mass and the momentum transfer squared of the isosinglet B20 GFF in 0{p^ 
is given by [65l | 



niN 



A. 



0,u+d lu+d 
20 "-B 



where B. 



0,u+d _ 
20 — 



B2o+'*(t = 0, = 0), and the terms AS; 



20 



dg t and 0^ 



CBChPT 



(28) 

are of 0{p^) and represent only 
, 0^ and dg may be obtained irom 
turns out to be unstable and produces large 
values for the counter term parameter ~ X5. This can be seen as indication that other higher order correction 

terms of 0{p^) not yet included in Eq. (pS)) are numerically important and needed to stabilize the extrapolation. We 
note that the counting scheme of ^65*1 suggests that Ai?2o is not a dominant C'(p^)-contribution concerning the pion 
mass dependence, at least for t — 0. This can be seen to some extent from the heavy-baryon-limit of Eq. (12811 . which 



a part the full 0{p ) contribution. The a priori unknown constants 

a fit to the lattice data. A fit to our lattice results based on Eq. ((28)) 



does not reproduce the full coefficient, cx (^4 + i3)2o'^'^, of the log(mJ)-term in HBChPT (see e.g. [H, [g^, 



but rather gives a term cx A2'^'^'^m'^log{rn^) without the -820"^'* 



m^log(m^) contribution. 
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FIG. 31: Lattice results for B^+'^ at t ^ -0.24 
GeV^ versus together with the result of a global 
simultaneous chiral fit based on Eqs. (I27|l , (I29|l and 
a variant of Eq. pSI) . as described in the text. 



FIG. 32; Lattice results for B^+''{t) at ^ 350 
MeV versus {—t) together with the result of a global 
simultaneous chiral fit based on Eqs. (|27p . (|29p and 
a variant of Eq. (I28|l . as described in the text. 



Since the instability of the fit can be traced back to the term Ai?2o^ (t,TO^), we performed the final fit dropping 
this contribution but keeping the counter terms cx t and cx m^. Based on this approach, we find that a global 
simuhaneous fit to the GFFs B^+'^ and 0^+^^, using Eqs. (I^T)) and (P^)) as described in the previous section, 

leads to a stable chiral extrapolation of all three GFFs. In particular the counter term parameters J^"^"* and 
gm„,u+d very small and fully compatible with zero within errors. We obtain i?2o"^'^ = —0.140(84) and 

^■zo'^'i^ ~ 0,TO^_phys) = —0.095(86) at the physical pion mass. Results of the fit are shown in Figs. (|3T|) and (|32|l . 
We cannot rule out that the inclusion of the full 0{p^) contributions to i32o^'* will lead to a qualitatively different 
dependence on t and in the region where lattice results are available, so that the results above should to be taken 
with due caution. 

In section |VH[ we will study the combined {t, m7r)-dependence of i?2o''' in HBChPT at 0{p^) [13]. One-loop graphs 
with insertions of pion operators are fully included, leading to a non-analytic dependence on the pion mass and the 
momentum transfer that is quite different from that shown in Figs. ([5T|) and ([5^ . 



F. CBChPT extrapolation of C^ifit) 
The (i,m^)-dependence of the isosinglet GFF C20 in CBChPT to 0{p'^) is given by ^ 

C-+\t, m^) = 1^Il^c%-+d + A'^^+\l+\t, m,) + AC2V'(i, m^) + 0{p^) , (29) 
niN 

where Cl^^"^ = €^^'^{1 = 0, = 0), and the term AC20+'' « i^Yu+d is a part of the fuU O(p3)-corrections [651]. 
In this counting scheme, counter terms of the form S^f^^'^ t and 6™""^^'^ rn^ first appear at 0{p'^). In order to get 
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a first idea about the possible t- and m^-dependence of in CBChPT, we have included the formally higher 

t,u+d_^ ;;nix,«+<i^2 ^^ic fit to our lattice data, resulting in a stable chiral extrapolation. 



order counter terms t and (5p 

With (x):; 



0.5 from Table Hill as an input parameter, and taking into account the CBChPT results for ^20''* ^^'^ 
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33: Lattice results for 
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^ versus together with the result of a global 
chiral fit using Eqs. J^HJ and 



FIG. 34: Lattice results for C^o'^it) at m„ ^ 350 
MeV versus (— t) together with the result of a global 
simultaneous chiral fit using Eqs. i27} . I|28p and 



^20'^ discussed in sections IV Dl and IV El respectively, we obtain from a global simultaneous fit C20"''' = —0.317(59) 
for the forward value in the chiral limit, and C2(^'^{t = 0, m7r,phys) = —0.267(62) at the physical pion mass. The 
parameters (5^"''''^ and turn out to be small. Changing (x)'^'^^ by ±10% results in a variation of C2i^'^'^{t = 0) 

by ±5%, which is significantly smaller than the statistical error of 19%. The corresponding extrapolations are shown 
in Figs. ([331) and ([Ml) . These results indicate a non-trivial dependence of the n — 2 moment of the isosinglet GPDs 
H and E on the longitudinal momentum transfer ^. 

A three-dimensional plot showing the combined (i, m7r)-dependence of Cjq^'' is presented in Fig. ([SSl) . where the 
error bars of the lattice data points are illustrated by the stretched cuboids. The lattice data are superimposed with 
the result from the chiral fit discussed above, which is shown as a surface. The statistical error bars of the fit are 
shown for clarity only as bands for t = and m^r = 0, respectively. In section fV J! below, we will compare the results 
based on CBChPT with a fit to C^r^"^ in the framework of HBChPT. 



G. Quark angular momentum J in CBChPT 

The forward limit values of the isovector and isosinglet GFFs A2o{t = 0) and B2o{t = 0) we have studied in 
sections IV Ai IV Bl IV Dl and IV El allow us to compute the angular momentum contributions of up- and down-quarks 
to the spin of the nucleon, J'' = l/2(yl2o(0) + i?|o(0)) — ^I'^ii^)'^ + -^2o('^))- From the separate chiral extrapolations 
of the isosinglet A20 and B20 in CBChPT, we find for the total u + d quark angular momentum at the physical 
pion mass J"^''(m7r,phys) = 0.213(44), corresponding to 43% of the total nucleon spin S = 1/2. Together with the 
chiral extrapolations for the isovector u — d combination, we obtain the surprising result that the total quark angular 
momentum is carried by the up-quarks, J"(m7r,phys) = 0.214(27) and that the contribution from down quarks is zero 
>/'^(^7T-7r,piiys) — —0.001(27). From Fig. [121 we note that the cancellation of and L"^ appears to be systematic 

for all m.^., and it will be interesting to understand whether this is accidental or has a physical origin. Taking into 
account preliminary results for the quark spin A^Q/2(t = 0) = AE^/2, as obtained from a ChPT extrapolation 
including the A resonance [l^llOl of gA = AE"^*^ and a self-consistently improved one-loop ChPT extrapolation of 
[4^^ -yirg find that the quark orbital angular momentum L'' = J'^ — AE''/2 contributions to the nucleon spin 
are L"^ = —0.195(44) and L"^ = 0.200(44) at the physical pion mass. The nearly complete cancellation of up and down 
quark 0AM contributions that we observe for pion masses above 350 MeV therefore also holds at at mjr^phys, where 
we find = 0.005(52). We emphasize again that no phenomenological values for AS = (1)aij, {x)q and {x)Aq 

have been included in the extrapolations, and that we have omitted disconnected diagrams in the lattice calculations. 
We will compare these CBChPT results with corresponding HBChPT results below in section FVJl 
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FIG. 35; Combined (i;,m^)-dependence of C2q'^ from a global simultaneous chiral fit using Eqs. \21\ . p8|l and (|29|) compared 
to lattice data. 



H. HBChPT extrapolation of and 

In heavy baryon chiral perturbation theory [g^, [g^] to 0{p'^), the combined (t, TO7r)-dependence of the GFF- 
combination E^^'^{t) = + t/ {^m\)Bl^+'^(t) is quite different from that of M^+'^{t) = A'^+'^{t) + B^+'^{t), 

which in the forward limit is equal to two times the total quark angular momentum 2Jq = M'^^'^it — 0). While at 
this order Ml^^'^ shows a non-analytic dependence on t and TOtt as discussed below, E'^q'^ is constant up to analytic 
tree-level contributions, 

El^+\t, m^) = i?°6"+' + EZ-^'-'i^l + E',t't. (30) 
A fit to our lattice results based on Eq. ([50]) is shown in Figs. ([55]) . (|57|) . ([55)1 and ([55]) . The linear dependence 




FIG. 36: Lattice results for at t = versus 

m% together with the result of a global chiral fit 
using Eq. dSOj. 



FIG. 37: Lattice results for £20^'* at -t « 0.24 
GeV^ versus together with the result of a global 
chiral fit using Eq. ((30)l . 



on t and works well even beyond the fitted range, i.e. for —t > 0.48 GeV^. This is not surprising since i?20^'' 
is clearly dominated by the GFF AI^q"^, which does not show a strong curvature in t as seen in Figs. ([8]) and pop . 
However, it is obvious that the HBChPT result in Eq. ([50)1 lacks structures which would allow for an upwards bending 
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of E'^+'^it = 0) = A'^+'^it = 0) at small pion masses towards the phenomenological value, in contrast to the covariant 
approach studied in section IV Al The fit gives 

^a«+<i ^ 0.481(15) in the chiral limit, which we will use for the 
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FIG. 38: Lattice results for E^^"^ at « 350 
MeV versus —t together with the result of a global 
chiral fit using Eq. (|30)l . 



FIG. 39: Lattice results for E^+'^ at m„ « 496 
MeV versus —t together with the result of a global 
chiral fit using Eq. ((30)l . 



chiral extrapolations based on HBChPT of the total angular momentum and the anomalous gravitomagnetic moment 
B^o'^it = 0) below. At the physical pion mass, we find i?20^''(t = 0,m^,phys) = {x)u+d = 0.485(14), which is 
approximately 10% below the phenomenological results, {x)2ld^'^ = 0.537(22) and {x)^^^™°^ = 0.538(22) [8l|. 
The pion mass dependence of Af^Q'"'^(t) for non-zero t is given by [6^. [6^ 



0,u+d 
20 



1 - 



(4^A)2 1^ A2 



(^) } + Aff ''^^(t + M^,-^^+%4 + M*-+^t, (31) 



The non-analytic dependence on t and m^r in M^''^\t,mTr) results 
from pion-operator insertions and is directly proportional to the (isosinglet) momentum fraction of quarks in the pion 



with new counter terms M. 
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and • 



in the chiral limit, ix) 



u+d 



We use (x] 



u+d 



0.5 from Table Hill for the fit. No additional parameters are needed to this 



order. The results of chiral fits based on Eq. ((3T|) are presented in Figs. (|40|) and (|4T 



We find 



0.522(41) 
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FIG. 40: Lattice results for M, 



u-\-d 



at \t\ 



0.24 



GeV^ versus together with the result of a global 
chiral fit using Eq. (|3TJ. 



FIG. 41: Lattice resuhs for M:^,^'' at 



350 



MeV versus {—t) together with the result of a 
global chiral fit using Eq. (|3ip. 



and M^+'^it = 0,TO^,phys) = 0.526(48). 



I. HBChPT extrapolation of B. 



u + d 



Total momentum and angular momentum conservation implies that the total anomalous gravitomagnetic moment 
of quarks and gluons in the nucleon has to vanish, ^„„B2o{t = 0) = 0. An interesting question is whether the 



individual quark and gluon contributions to i?2o are separately zero or very small, as previously speculated [82 
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The first lattice calculations [T3, [3 showed that BI^q'^ is compatible with zero for (— t) > 0.5 GeV^ at pion masses 



of TOtt > 600 MeV. Based on our new results and the ChPT fits performed above, we are now in a position to study 
^20'^ more carefully as a function of t and m^r . The GFF -620^'' can be written as a linear combination of Eq. ([50)) and 
Eq. dnU. A separate fit to the data with fixed ^^20'^'' = 0.481(15) gives B^+'^{t = 0,TO^,phys) = 0.050(49), which is 
compatible with the fits to M^^"^ and £'20^'' above that in combination give {M - E)1^'^(t = 0, m^^phys) = 0.041(50). 
Although the absolute value of B'l^^^it = 0) is again rather small, we note that the sign is different from that found 
in section FV El based on the CBChPT fit. A 10% variation of the input parameter {xfj^^ leads to change of 0.023 in 

^20^'*(^i 'Tijr.phys) at i = 0, which is well below the statistical error of 0.049, and a change of 0.008 at a momentum 
transfer of \t\ k, 0.24 GeV^, which is well below the statistical error of 0.031. 

The dependence of B'^q'^ on t and on the pion mass is shown in Figs. and The dependence on the momen- 

) where contributions 

have not been included, but the two results are statistically consistent. The 



tum transfer squared turns out to be somewhat different from the CBChPT result in Fig. 
from pion operator insertions oc 
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FIG. 42: Lattice results for Bgo^'' at \t\ « 0.24 
GeV'^ versus together with the result of a global 
chiral fit using a linear combination of Eqs. (|30|l and 
(ED- 



FIG. 43; Lattice results for Bjo^'' at Ri350 MeV 
versus (— together with the result of a global 
chiral fit using a linear combination of Eqs. pO|l 
and (PT|) . 



combined dependence of B^^''' on i and m-p; from the HBChPT fit compared to lattice data points is presented in 
a 3d-plot in Fig. (|^l)l . It is interesting to note that the fit based on Eq. (|3ip leads to a clearly visible non-analytic 
dependence of B^^'^ on the pion mass and the momentum transfer. In particular, wc find a non-zero, neg ative Bl^'^ 
for \t\ > 0.05GeV^, ^ O.lGeV'^ from the chiral extrapolation, which may be confirmed in future lattice calculations 
or experimental measurements. 



J. HBChPT extrapolation of C^^"^ 

At order 0(jp-), the pion mass dependence of the GFF C^^'^{t) is given by [H,!!!,!!!] 

Cl^\t,m^) = 1 _ {^"""^' + """^ + ""-^ + C""-" +'^m^ + C*'"+'^ ^| , (32) 

where Cjq"^'' = C^^'^it — 0,70^ — 0). The terms i?2^''^''(<:, m^) and E^r^'^\t,mTr) contain non-analytic terms in t and 

that come from insertions of pion operators proportional to (x)'^'^^. Additionally, E^''^^ {t, mT^) depends on the low 
energy constants ci, C2 and C3. We fix these parameters to the values given in Table Hill The result of a fit to our lattice 
data as a function of the pion mass squared for fixed t ~ —0.24 GeV^ is shown in Fig. dS) The t-dependence at a pion 
mass of w 350 MeV is presented in Fig. (06]). We find C",^ ^'^ = -0.507(55), and Cl^^^'^it = 0, m^,phys) = -0.421(54) at 
the physical pion mass. We note that these values are approximately 60% larger in magnitude than the corresponding 
CBChPT results in section IVFl based on covariant ChPT, which is directly related to the stronger downwards bending 
of Cjo^'' for {-t) ^ in Fig. SS] compared to the slight upwards bending in Fig. (|33p . As in the case of the CBChPT 
extrapolation, a variation of the input value by 10% results in a w 5% change of C2o"^'^(t — 0), which is below 

the 11% statistical error. 

Fig. (|47)) shows our combined lattice results for C20^'^ versus (— i) and TOtt in a single three-dimensional plot, together 
with the result from the HBChPT fit discussed above, which is shown as a surface. The statistical error bars are 
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FIG. 44: Combined {t, m.^)-dependence of the quark anomalous gravitomagnetic moment B2q'^ from a global chiral fit using a 
linear combination of Eqs. (I30|l and (I3ip compared to lattice data. 
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FIG. 45: Lattice results for C. 



at t\ 



0.24 



GeV^ versus together with the result of a global 
chiral fit using Eq. p2)l. 



FIG. 46; Lattice results for C^^"^ at m^^350 MeV 
versus — t together with the result of a global chiral 
fit using Eq. (|32)l . 



shown for clarity as bands for t — GeV^ and = MeV, respectively. It is interesting to note that the overall 
shape of the extrapolation surface is similar to the CBChPT result in section IV Fl The only behavior that differs 
by more than the statistical errors is the slightly stronger bending towards negative values of C^q'^ at the origin in 

Fig. dlZD. 



K. HBChPT extrapolation of quark angular momentum J 



From our results for M20 above, we find a total quark angular momentum J"+'^ = 0.263(24) at the physical pion 
mass, which is larger than but statistically compatible with the CBChPT value in section fV Gl 

As an alternative, we can also calculate J'^+'^ by first extrapolating B2o{t^m-^) to t = and combining it with 
^20 = 0, vItt) to obtain J(m.n.), and then extrapolating the values of J(m,r) to the physical pion mass using HBChPT 
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FIG. 47: Combined (t, m^)-dependence of Cjo from a global chiral fit using Eq. (l32ll compared to lattice data. 



that explicitly includes the A resonance [55] . Evaluating Eq. PT|) at t = Q yields 



1 



J' 



(33) 



which agrees with |i55j . Note that in the notation of [55| we have bgN = (A + 5)20"+ = mO^"+'^. Including explicitly 



0.3 

0.25 
0.2 

0.15 
0.1 

0.05 



0.1 



0.2 0.3 0.4 
[GeV^] 



0.5 0.6 



FIG. 48: Chiral extrapolation of J""'"'' usmg 
HBChPT including the A resonance, Eq. (f5I)) . 
The fit and error band on the axis are explained 
in the text. 
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FIG. 49: Comparison of a global fit to ^30+ +^30+ 
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fit and error band on the axis are explained in the 
text. 



the A resonance in the calculation, the ChPT result then reads [5£ 



r+'im^; A) = r+'^M - i Q(A + B) 
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(34) 



29 



where A — mi^ — mjsi denotes the A-nucleon mass difference. In order to reduce the number of free parameters 
in the fit, we use A = 0.3 GeV and the large- iVc relation gT^NA = 3/(2'^/^)5^ from Tabic Hill We first extrapolate 
B^o'^it^mj,) linearly in i to t = 0. A fit based on Eq. ^ to the lattice results for {A + B)'^^'^{t = 0, m^) with 
< 700 MeV then gives {A + Bf^^+'^ = bqN = 0.545(12), 6,a = 0.427(51) and J"+'^(m^,phys; A) = 0.212(32) 
at the physical pion mass, which is very close to the CBChPT result in section IV Gi J"^'^(m^,phys) = 0.213(44). 
The result for J"+'^(m^; A) as a function of the pion mass is shown in Fig. (pS)) . where the error due to the linear 
extrapolation of B{t,mTr) to t = is indicated by the error band at the mj-axis as explained at the end of section 
IIVI and the error bars on the lattice data points for J only include the errors arising from A2o{t = 0). In Fig. (j49p . 
we compare the result of the chiral extrapolation of M2f^'^(t) for non-zero t from the previous section with the lattice 
results for J^+^ corresponding to the extrapolated B^^'^lt ~ 0). The two different ways of fitting and extrapolating 
J"+'^ in are compatible within errors, Jf^tmt^o = 0.263(24) versus J^J^g (m^,phys; A) = 0.212(32), where the 
chiral fit including the A leads to a stronger curvature at small and therefore to a smaller central value for 
J"+'^ at the physical point. Together with prehminary results for the quark spin A"Q'^/2{t = 0) = AE"+'^/2 as 
obtained from a self-consistently improved one-loop ChPT extrapolation [48|, we find that the quark orbital angular 
momentum L'^+d- = j"-i-<i _ AE"+''/2 contribution to the nucleon spin is zero within errors: = 0.005(43) for 

Jt^o{m^,phys;A) = 0.212(32) and = 0.056(37) for J^+^ = 0.263(24) at the physical pion mass. 



VI. SUMMARY AND CONCLUSIONS 



This work presents the first comprehensive full lattice QCD study of the lowest three moments of unpolarized and 
polarized GPDs in the chiral regime with pion masses as low as 350 MeV. We find good overall agreement with existing 
experimental results. We note, however, that we have omitted disconnected diagrams, which in principle contribute 
to isoscalar observables. 

As in our previous study of the axial vector coupling constant [85], the consistency of these moments in lattice 
volumes of (2.5fm)'^ and (3.5fm)'^ at = 350 MeV indicates that finite volume effects are negligible within statistical 
errors. 

One significant result of this work is the clear indication that the transverse size of the nucleon, as characterized 
by the transverse 2d rms radius, is a strongly decreasing function of the momentum fraction x carried by the quarks. 
At the lightest quark mass, the isovector transverse rms radius drops by almost 60% between the zeroth moment, 
which roughly corresponds to an average momentum fraction[86] (x) « 0.2 and the second moment, which roughly 
corresponds to an average momentum fraction (x) « 0.4. This decrease in the chiral regime is even stronger than our 
original observation of the decrease of the transverse size with momentum fraction in the "heavy pion world" [U ■ 

In a first direct comparison with phenomenological parametrizations of the GPDs H{x, £, = Q,t) and H{x, £^ — Q,t) 
constrained by structure function and form factor data ll| , we find qualitative consistency for the ratios of GFFs in 
both the isovector and isosinglet cases. 

Our results provide insight into the contributions of the spin and orbital angular momentum of u and d quarks to 
the spin of the proton. Although the individual orbital angular momentum contributions of the u and d quarks are 
sizeable, L"^ w — « 30%, they cancel within errors so that the total contribution is « 0. In addition, the spin 

and orbital contributions of the d quark also cancel within errors, so J"* « 0. The total quark angular momentum 
contribution is J"+'^ « 40 — 50% at our lowest pion masses. 
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J" 




covariant BChPT 

HBChPT 
HBChPT with A 
phenomenology 


0.157(10) 
0.155(5) 


0.273(63) 


-0.017(41) 


0.520(24) 
0.485(14) 

0.537(22) 


-0.095(86) 
0.050(49) 


-0.267(62) 
-0.421(54) 


0.213(44) 
0.263(24) 
0.212(32) 


0.214(27) 


-0.001(27) 



TABLE V: Summary of proton observables &t t = Q and phys from chiral extrapolations, in the MS scheme at a scale of 



More quantitatively, we performed a variety of chiral fits to the unpolarized n — 2 moments using covariant BChPT 
[65l | and HBChPT results [H, [g^, [gl, [g^l • A summary of our results for various observables at the physical pion mass 
and vanishing momentum transfer t = is given in Table |Vl where the quoted errors are statistical only. We note that 
a consistent inclusion of all 0{p'^) terms in our fits, which may involve additional constants, would have the potential 
to increase the statistical errors on the physical quantities at the chiral limit. However, for the reasons described in 
section IV A[ the prescription of adding a single m\ term with extremal values of "natural" coefficients does not lead 
to reasonable fits, so we have not attempted to include quantitative estimates of these errors in the present work. 
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With the exception of a fit to Jq based on HBChPT including the A resonance, we have consistently extrapolated 
the GFFs simultaneously in the pion mass and the momentum transfer squared t. The simultaneous covariant BChPT 
fits to the GFFs A2Q, B20 and C20, which include approximately 120 lattice data points and typically 9 unknown 
low energy constants, produce reasonable parametrizations of the {t, m^)-dependences of the generalized form factors 
in the ranges < 700MeV and \t\ < 0.5 GeV^. In particular, the covariant extrapolations for the isovector 
and isosinglet momentum fractions (x) yield values at the physical point remarkably close to phenomenology. This 
represents a significant advance in our understanding of the pion mass dependence of these important observables. 
The first exploration of the combined non-analytic dependence of the isosinglet GFFs -620^'* and Cjq'''^ on t and TO,r 
was made using covariant BChPT and HBChPT, and visualizations of the resulting (i, m7r)-dependences of these 
GFFs in 3d-plots reveal interesting non-linear correlations in the pion mass and the momentum transfer squared. 

In spite of the overall success of the chiral extrapolations, it is clear that the ChPT analysis has not yet been 
carried to sufficiently high order to be applicable to the full range of pion masses included in this work. The facts that 
HBChPT fits to A^q"^ and A^q"^ cannot describe the behavior of the lattice data over as large a range or as accurately 
as covariant BChPT, and that the fitted counter terms are so different indicate that the higher order terms in l/rriN 
included in CBChPT are important for these observables. Similarly, the fact that it was important to include some 
particular terms of order 0{p^) and 0{p'^) in some CBChPT fits indicates the need to fully determine these orders in 
ChPT so that they can be consistently included in fits to lattice data. Finally, the significant effect of including the A 
in J"+'*, its known importance in the axial charge, and large Nc arguments [8 Tj indicate the desirability of consistent 
inclusion of the A. Thus, future progress requires both the extension of lattice calculations to lower pion mass and 
the inclusion of higher order terms and A degrees of freedom in ChPT. 
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APPENDIX A: TABLES 



Lattice parameters of the datasets 
scale of /x^ = 4 GeV^. 



1 , . . . , 6 are provided in Table [J The GFFs are given in the MS scheme at a 



-t[GeV^ 



AU-, 

^10 



r)U — 



AU-, 

^20 



-"20 



Cu — 
20 



^30 



ryu — 
^30 



0.000 
0.203 
0.243 
0.243 
0.420 
0.475 
0.477 
0.498 
0.697 
0.741 
0.911 
0.918 
0.996 
1.117 
1.199 
1.239 



1.000(1) 

0.731(73) 

0.796(6) 

0.735(29) 

0.606(64) 

0.659(9) 

0.605(29) 

0.646(15) 

0.561(12) 

0.542(16) 

0.487(17) 

0.462(31) 

0.493(16) 

0.425(17) 

0.412(21) 

0.426(15) 



2.331(347) 

2.645(55) 

2.570(120) 

1.748(239) 

2.093(48) 

2.000(110) 

1.994(60) 

1.665(51) 

1.611(57) 

1.434(65) 

1.317(125) 

1.289(47) 

1.152(58) 

1.100(70) 

1.077(42) 



0.223(3) 

0.179(17) 

0.201(3) 

0.187(7) 

0.160(16) 

0.180(3) 

0.170(7) 

0.179(4) 

0.163(4) 

0.160(5) 

0.149(6) 

0.136(9) 

0.154(5) 

0.137(6) 

0.129(6) 

0.137(5) 



0.331(91 
0.398(16 
0.340(29 
0.258(53 
0.341(15 
0.311(29 
0.325(16 
0.294(15 
0.288(16 
0.292(23 
0.226(33 
0.254(16 
0.233(18 
0.225(21 
0.214(13 



-0.094(113) 

-0.002(11) 

-0.006(28) 

-0.065(32) 

-0.008(6) 

-0.027(14) 

0.014(7) 

-0.001(6) 

0.006(6) 

0.005(8) 

0.012(15) 

0.006(8) 

-0.001(5) 

0.014(7) 

0.004(5) 



0.072(4) 

0.072(10) 

0.071(2 

0.066(3 

0.055(7 

0.066(3 

0.062(3 

0.066(3 

0.061(3 

0.059(3 

0.055(4 

0.053(4 

0.053(4 
0.047(4 
0.054(4 



0.051(59) 

0.109(9) 

0.078(13) 

0.078(29) 

0.092(8) 

0.069(13) 

0.088(9) 

0.084(8) 

0.075(8) 

0.102(12) 

0.057(15) 

0.070(5) 

0.074(9) 

0.073(10) 

0.069(7) 



TABLE VL Results for the isovector unpolarized generalized form factors for dataset 1. 



-t[GeV'l A'lo 



10 



0.000 
0.200 
0.242 
0.243 
0.414 
0.473 
0.475 
0.498 
0.693 
0.741 
0.904 
0.911 
0.996 
1.106 
1.195 
1.238 



1.000(2) 

0.924(100) 

0.798(8) 

0.804(38) 

0.782(85) 

0.678(11) 

0.656(37) 

0.651(17) 

0.587(16) 

0.549(18) 

0.506(23) 

0.556(56) 

0.464(20) 

0.450(23) 

0.437(29) 

0.404(18) 



3.454(490) 

2.540(67) 

2.384(146) 

2.226(301) 

2.034(57) 

1.697(133) 

1.873(65) 

1.682(63) 

1.441(57) 

1.391(80) 

1.361(202) 

1.189(51) 

1.147(68) 

1.085(88) 

0.987(41) 



Aii.-d 

^20 

0.223(4) 

0.196(26) 

0.201(4) 

0.196(9) 

0.189(22) 

0.181(4) 

0.171(9) 

0.177(5) 

0.166(6) 

0.153(5) 

0.150(8) 

0.160(15) 

0.146(7) 

0.137(8) 

0.143(9) 

0.134(6) 



riu — 

-^20 



Cu — 
20 



0.461(125) 

0.350(20) 

0.289(36) 

0.346(69) 

0.314(18) 

0.221(35) 

0.274(18) 

0.297(20) 

0.222(17) 

0.224(24) 

0.218(49) 

0.212(18) 

0.210(20) 

0.226(27) 

0.189(14) 



0.261(149) 

-0.006(14) 

-0.044(37) 

0.005(43) 

-0.001(7) 

0.000(19) 

-0.008(10) 

0.001(8) 

-0.007(8) 

-0.002(11) 

0.007(22) 

-0.005(10) 

0.005(7) 

-0.010(11) 

-0.001(6) 



AU-d 
^30 

0.078(6) 

0.054(15) 

0.070(3) 

0.072(5) 

0.061(10) 

0.063(4) 

0.066(5) 

0.064(4) 

0.062(5) 

0.059(4) 

0.058(6) 

0.066(8) 

0.049(5) 
0.052(5) 
0.055(5) 



-^30 



0.168(76) 

0.072(11) 

0.057(17) 

0.137(34) 

0.073(10) 

0.053(18) 

0.050(10) 

0.077(11) 

0.045(9) 

0.067(15) 

0.039(23) 

0.058(6) 

0.052(11) 

0.044(12) 

0.045(9) 



TABLE VII: Results for the isovector unpolarized generalized form factors for dataset 2. 



32 



-t[GeV^' 



0.000 
0.193 
0.241 
0.242 
0.403 
0.469 
0.472 
0.498 
0.685 
0.740 
0.891 
0.901 
0.996 
1.089 
1.189 
1.238 



AU — d 

1.000(1) 
0.674(83) 

0.770(8) 

0.724(35) 

0.627(78) 

0.640(12) 

0.584(34) 

0.593(17) 

0.546(17) 

0.495(20) 

0.439(20) 

0.397(41) 

0.426(19) 

0.391(23) 

0.359(26) 

0.356(18) 



TJU — 

-°io 



2.147(465) 

2.399(62) 

2.232(154) 

1.723(305) 

1.890(54) 

1.713(142) 

1.819(61) 

1.543(62) 

1.412(65) 

1.101(66) 

1.065(180) 

1.131(51) 

0.938(63) 

0.867(71) 

0.939(43) 



^20 

0.226(4) 

0.176(24) 

0.200(4) 

0.187(8) 

0.165(20) 

0.181(4) 

0.171(9) 

0.169(5) 

0.165(6) 

0.152(6) 

0.136(7) 

0.128(13) 

0.142(7) 

0.124(8) 

0.127(9) 

0.125(7) 



JDU — 

-°20 



0.226(119) 

0.358(18) 

0.345(35) 

0.243(68) 

0.300(17) 

0.265(34) 

0.305(19) 

0.263(20) 

0.249(19) 

0.186(20) 

0.155(42) 

0.232(17) 

0.172(17) 

0.143(22) 

0.207(14) 



'^20 



-0.044(152) 

-0.009(13) 

0.052(35) 

-0.010(45) 

0.000(7) 

0.003(20) 

0.013(9) 

-0.007(9) 

0.007(9) 

0.015(11) 

0.034(19) 

-0.001(9) 

0.011(7) 

-0.001(10) 

-0.004(6) 



A u — d 

0.074(5) 

0.046(13) 

0.074(3) 

0.067(4) 

0.049(9) 

0.069(4) 

0.065(5) 

0.064(4) 

0.065(5) 

0.058(4) 

0.045(5) 

0.047(6) 

0.053(5 

0.049(5 
0.047(5 



30 



0.058(87) 

0.094(12) 

0.075(17) 

0.118(38) 

0.078(10) 

0.075(17) 

0.087(11) 

0.075(11) 

0.075(10) 

0.039(15) 

0.044(22) 

0.053(7) 

0.046(12) 

0.046(11) 

0.055(9) 



TABLE VIII: Results for the isovector unpolaxized generalized form factors for dataset 3. 



-t[GeV" 



"^10 



Bo, 



^3(1 



0.000 
0.188 
0.240 
0.241 
0.394 
0.465 
0.469 
0.498 
0.678 
0.739 
0.880 
0.892 
0.996 
1.072 
1.183 
1.237 



1.000(2) 

0.548(129) 

0.753(13) 

0.765(60) 

0.529(110) 

0.601(19) 

0.630(57) 

0.595(28) 

0.505(27) 

0.501(30) 

0.420(38) 

0.549(115) 

0.400(28) 

0.376(34) 

0.462(58) 

0.347(25) 



3.044(876) 

2.281(96) 

2.335(250) 

1.927(522) 

1.696(83) 

2.034(223) 

1.642(99) 

1.306(95) 

1.273(94) 

1.113(126) 

1.470(396) 

1.018(70) 

0.840(96) 

0.800(143) 

0.737(54) 



0.214(7) 

0.202(46) 

0.190(6) 

0.182(14) 

0.143(31) 

0.173(7) 

0.161(14) 

0.171(8) 

0.159(9) 

0.147(8) 

0.138(14) 

0.152(32) 

0.124(10) 

0.115(12) 

0.140(17) 

0.118(9) 



0.505(205) 

0.307(30) 

0.309(56) 

0.188(118) 

0.271(26) 

0.311(52) 

0.272(28) 

0.208(29) 

0.224(28) 

0.196(37) 

0.202(83) 

0.159(26) 

0.148(30) 

0.159(43) 

0.163(21) 



-0.059(262) 

-0.015(21) 

0.091(57) 

0.119(81) 

-0.020(12) 

0.003(32) 

-0.011(15) 

-0.024(13) 

0.015(14) 

0.012(17) 

-0.036(45) 

0.007(13) 

0.010(11) 

0.000(17) 

0.007(9) 



0.070(8) 

0.074(28) 

0.071(5) 

0.066(6) 

0.047(16) 

0.070(6) 

0.066(7) 

0.067(6) 

0.063(8) 

0.061(6) 

0.062(11) 

0.062(15) 

0.055(8) 

0.057(10) 

0.041(8) 



0.045(165) 

0.059(20) 

0.095(28) 

0.001(71) 

0.072(16) 

0.083(28) 

0.093(15) 

0.072(19) 

0.063(14) 

0.032(28) 

0.052(39) 

0.030(11) 

0.052(21) 

0.066(22) 

0.048(15) 



TABLE IX: Results for the isovector unpolarized generalized form factors for dataset 4. 



u—<t 
9n 



-t[GeV' 



^10 



ryu — 
-°10 



^20 



TJU — 

-"20 



AU-. 

^30 



30 



0.000 
0.177 
0.238 
0.240 
0.378 
0.459 
0.464 
0.498 
0.665 
0.738 
0.858 
0.875 
0.996 
1.042 
1.173 
1.235 



1.000(3) 

0.387(251) 

0.723(24) 

0.587(97) 

0.604(35) 

0.502(100) 

0.466(42) 

0.491(51) 

0.368(54) 

0.349(67) 



0.205(14) - 

2.385(1.795) 0.095(79) 0.717(475) 0.129(518) 



1.957(166) 
1.541(396) 

1.566(137) 
1.072(373) 
1.191(149) 
1.228(166) 
1.091(156) 
0.540(208) 



0.172(11) 0.282(49) 
0.135(20) 0.262(91) 



0.156(12) 
0.131(25) 
0.127(12) 
0.134(17) 
0.126(15) 
0.115(25) 



0.221(44) 
0.344(97) 
0.152(40) 
0.219(51) 
0.206(46) 
0.203(74) 



-0.025(36) 
-0.095(93) 

0.009(22) 

0.047(65) 

0.001(29) 

-0.034(27) 

0.071(32) 

0.010(34) 



0.063(16) - 

-0.033(63) 0.537(408) 
0.049(11) 0.020(40) 
0.063(12) 0.045(52) 



0.043(12) 
0.052(15) 
0.037(10) 
0.004(18) 
0.033(11) 
0.012(24) 



0.377(60) 0.606(121) 0.135(22) 0.166(50) -0.009(25) 



0.270(75) 0.411(176) 0.112(25) 0.057(53) -0.054(29) 0.036(13) 
0.291(42) 0.589(88) 0.090(16) 0.074(31) -0.019(14) 0.027(16) 



0.065(33) 
0.069(52) 
0.025(29) 
0.011(40) 
0.045(30) 
0.048(61) 

0.036(22) 

0.039(34) 
0.014(25) 



TABLE X: Results for the isovector unpolarized generalized form factors for dataset 5. 



33 



-t[GeV" 



^1 



TJU — 

-°io 



A-u-d 
^20 



20 



20 



^30 



T>U — 

^30 



0.000 
0.107 

0.124 
0.125 
0.219 
0.244 
0.245 
0.254 
0.359 
0.379 
0.471 
0.473 
0.508 
0.578 
0.615 
0.632 



1.000(4) 
1.035(192) 

0.850(19) 

0.829(84) 

0.868(190) 

0.767(28) 

0.733(79) 

0.678(39) 

0.693(36) 

0.584(40) 

0.617(42) 

0.551(91) 

0.602(56) 

0.561(45) 

0.478(49) 

0.512(45) 



3.055(997 
2.756(233 
2.262(376 
2.749(823 
2.330(193 
1.883(328^ 
2.074(219 
1.910(170 
1.786(192 
1.625(185 
1.910(448 
1.533(172 
1.447(150 
1.324(186 
1.268(131 



0.206(14 
0.190(45 
0.197(11 
0.192(20 
0.167(42 
0.177(11 
0.179(20 
0.179(14 
0.161(12 
0.158(14 
0.162(15 
0.155(28 
0.185(21 
0.155(15 
0.135(17 
0.173(18 



0.032(253) 

0.331(59) 

0.284(110) 

0.182(169) 

0.298(49) 

0.137(93) 

0.273(61) 

0.251(50) 

0.250(52) 

0.140(57) 

0.116(107) 

0.227(64) 

0.188(43) 

0.117(59) 

0.193(46) 



-0.083(435) 

-0.097(54) 

-0.017(113) 

0.026(130) 

-0.015(28) 

-0.164(73) 

0.004(43) 

0.005(33) 

0.010(39) 

-0.009(35) 

-0.030(56) 

-0.024(46) 

-0.028(20) 

-0.060(34) 

-0.040(25) 



0.078(16) 
0.060(32) 
0.070(10) 
0.072(12) 
0.063(24) 
0.063(10) 
0.068(12) 
0.081(13) 
0.061(12) 
0.067(11) 
0.064(14) 
0.077(17) 

0.068(13) 
0.068(13) 
0.076(17) 



-0.303(221) 

0.099(46) 

0.094(62) 

-0.084(109) 

0.091(34) 

0.099(64) 

0.067(42) 

0.087(36) 

0.040(36) 

0.019(50) 

0.065(67) 

0.044(26) 

0.031(34) 

0.012(37) 

0.002(36) 



TABLE XI: Results for the isovector unpolarized generalized form factors for dataset 6. 



u + d 



rju + d 
_°10__ 



u + d 



^20 



u+cL 
20 



^30 



-t[GeV^] A 



Al 



0.000 
0.203 
0.243 
0.243 
0.420 
0.475 
0.477 
0.498 
0.697 
0.741 
0.911 
0.918 
0.996 
1.117 
1.199 
1.239 



3.000(2) 

2.068(180) 

2.229(14) 

2.067(71) 

1.649(153) 

1.756(17) 

1.633(63) 

1.689(28) 

1.430(24) 

1.342(28) 

1.177(34) 

1.105(61) 

1.157(25) 

0.992(35) 

0.911(35) 

0.955(24) 



-0.828(581) 

0.067(67) 

0.083(150) 

-0.380(350) 

0.008(60) 

0.117(135) 

0.148(78) 

-0.026(64) 

0.147(72) 

0.049(78) 

0.112(156) 

0.098(64) 

0.024(71) 

0.086(80) 

0.090(52) 



0.590(5) 

0.473(40) 

0.506(5) 

0.473(15) 

0.398(35) 

0.437(5) 

0.415(15) 

0.428(7) 

0.383(7) 

0.371(8) 

0.345(9) 

0.316(17) 

0.343(8) 

0.303(10) 

0.281(10) 

0.301(8) 



-0.230(145) 

-0.024(20) 

-0.067(44) 

-0.178(91) 

-0.042(19) 

-0.044(43) 

-0.001(23) 

-0.048(20) 

0.000(21) 

-0.040(24) 

-0.067(44) 

0.003(24) 

-0.040(21) 

-0.016(24) 

-0.009(18) 



-0.201(196) 

-0.130(19 

-0.086(51 

-0.056(53 

-0.090(10 

-0.116(23 

-0.062(12 

-0.072(10 

-0.057(10 

-0.086(13^ 

-0.100(26 

-0.074(11 

-0.066(9) 

-0.053(11) 

-0.058(7) 



0.176(7) 

0.178(19) 

0.163(4) 

0.158(6) 

0.137(14) 

0.144(5) 

0.143(7) 

0.148(5) 

0.132(7) 

0.129(5) 

0.120(8) 

0.117(8) 

0.111(7) 
0.100(6) 
0.110(6) 



-0.123(95) 

-0.009(13) 

-0.024(20) 

-0.082(45) 

-0.013(12) 

-0.024(19) 

-0.007(14) 

-0.015(13) 

-0.009(12) 

0.002(16) 

-0.015(21) 

0.061(8) 

-0.009(13) 

0.009(13) 

-0.004(12) 



TABLE XIL Results for the isosinglet unpolarized generalized form factors for dataset 1. 



AU+d 
^10 



10 



\u+d 



■yu+d 



u + d 

on 



\u+d 



u+d 
30 



-t[GeV' 



0.000 
0.200 
0.242 
0.243 
0.414 
0.473 
0.475 
0.498 
0.693 
0.741 
0.904 
0.911 
0.996 
1.106 
1.195 
1.238 



3.000(3) 

2.286(222) 

2.227(17) 

2.134(91) 

1.773(178) 

1.763(23) 

1.666(86) 

1.684(36) 

1.430(34) 

1.336(38) 

1.229(48) 

1.193(119) 

1.094(33) 

1.023(49) 

0.925(57) 

0.916(31) 



0.641(651) 

0.049(88) 

0.186(183) 

0.352(392) 

0.055(75) 

0.191(179) 

0.001(98) 

0.073(79) 

0.049(88) 

0.009(99) 

-0.001(229) 

-0.023(73) 

-0.008(83) 

-0.067(120) 

-0.045(65) 



0.577(7) 

0.478(51) 

0.496(6) 

0.474(19) 

0.416(42) 

0.431(7) 

0.403(19) 

0.419(10) 

0.383(10) 

0.357(10) 

0.334(14) 

0.347(32) 

0.326(11) 

0.299(15) 

0.292(16) 

0.292(10) 



-0.010(188) 

-0.035(27) 

-0.009(50) 

0.083(95) 

0.000(23) 

-0.048(50) 

-0.050(27) 

0.013(25) 

-0.035(26) 

-0.026(31) 

-0.040(68) 

-0.017(27) 

-0.007(24) 

-0.021(36) 

-0.025(21) 



0.077(243) 
-0.119(25) 
-0.125(63) 
-0.087(67) 
-0.081(11) 
-0.093(30) 
-0.069(16) 
-0.068(12) 
-0.069(14) 
-0.045(18) 
-0.035(39) 
-0.061(13) 
-0.040(13) 
-0.033(17) 
-0.050(9) 



0.170(9) 

0.131(26) 

0.152(5) 

0.153(8) 

0.124(16) 

0.137(7) 

0.138(8) 

0.136(6) 

0.128(9) 

0.123(6) 

0.114(9) 

0.129(13) 

0.102(10) 

0.107(9) 

0.107(7) 



-0.029(131) 

-0.033(17) 

-0.022(25) 

0.064(49) 

-0.006(15) 

-0.014(26) 

-0.055(19) 

0.010(17) 

-0.040(16) 

-0.020(23) 

-0.046(35) 

0.059(10) 

-0.001(17) 

-0.031(20) 

-0.018(15) 



TABLE XIIL Results for the isosinglet unpolarized generalized form factors for dataset 2. 



34 



-t[GeV^ 



u-\-d 



ju-\-cl 



u-\-cl 
20 



^30 



■u+d 



A 



AU + . 

^20 



0.000 
0.193 
0.241 
0.242 
0.403 
0.469 
0.472 
0.498 
0.685 
0.740 
0.891 
0.901 
0.996 
1.089 
1.189 
1.238 



3.000(3) 

2.036(204) 

2.172(16) 

2.047(87) 

1.578(169) 

1.688(23) 

1.608(80) 

1.615(33) 

1.365(35) 

1.268(35) 

1.067(39) 

0.979(81) 

1.036(31) 

0.882(46) 

0.874(47) 

0.841(29) 



0.517(614) 

-0.034(82) 

0.024(170) 

0.071(377) 

-0.023(70) 

-0.156(161) 

0.002(81) 

-0.065(75) 

-0.048(81) 

-0.116(85) 

0.123(180) 

-0.024(67) 

-0.111(74) 

-0.098(91) 

-0.041(56) 



0.565(7) 

0.429(47) 

0.482(6) 

0.456(18) 

0.382(41) 

0.420(7) 

0.401(19) 

0.402(9) 

0.369(10) 

0.351(10) 

0.316(12) 

0.292(23) 

0.317(11) 

0.286(14) 

0.279(15) 

0.273(10) 



-0.111(175) 

-0.002(24) 

0.012(50) 

-0.040(95) 

-0.019(22) 

0.005(50) 

0.017(25) 

-0.010(25) 

-0.010(25) 

-0.057(28) 

-0.051(60) 

-0.011(25) 

-0.039(25) 

-0.033(30) 

-0.008(19) 



-0.034(247) 

-0.110(20 

-0.094(57; 

-0.113(71 

-0.092(11 

-0.049(29' 

-0.048(14' 

-0.072(13; 

-0.042(13' 

-0.054(16; 

-0.013(30' 

-0.050(12; 

-0.046(11 

-0.041(15 

-0.042(8) 



0.171(9) 

0.111(23) 

0.159(5) 

0.152(7) 

0.123(17) 

0.143(6) 

0.143(9) 

0.138(5) 

0.133(9) 

0.125(6) 

0.098(10) 

0.102(11) 

0.106(10) 

0.096(8) 
0.103(7) 



-0.086(131) 

-0.014(18) 

-0.012(25) 

0.001(54) 

-0.014(14) 

-0.014(26) 

0.003(16) 

-0.013(16) 

0.006(15) 

-0.052(21) 

0.007(32) 

0.056(10) 

-0.029(17) 

-0.003(17) 

-0.001(14) 



TABLE XIV: Results for the isosinglet unpolarized generalized form factors for dataset 3. 



^10 



^20 



20 



rtu + cL 
^30 



-t[GeV^ 



4^ 



0.000 
0.188 
0.240 
0.241 
0.394 
0.465 
0.469 
0.498 
0.678 
0.739 
0.880 
0.892 
0.996 
1.072 
1.183 
1.237 



3.000(3) 

1.817(355) 

2.127(26) 

2.098(153) 

1.382(262) 

1.632(35) 

1.611(130) 

1.521(55) 

1.317(50) 

1.225(59) 

1.090(79) 

1.246(232) 

0.960(45) 

0.878(68) 

0.960(113) 

0.789(39) 



-1.870(1.160) 

-0.167(128) 

0.165(296) 

-0.670(677) 

-0.085(113) 

-0.008(263) 

0.016(131) 

-0.180(116) 

0.033(131) 

0.081(152) 

0.252(387) 

0.026(101) 

-0.009(122) 

0.032(182) 

-0.024(81) 



0.530(12 

0.429(84 

0.459(9) 

0.442(30 

0.319(59^ 

0.398(10 

0.383(30 

0.385(15 

0.346(14 

0.334(15 

0.316(23 

0.375(68 

0.284(16 

0.245(20 

0.292(31 

0.246(14 



-0.034(281) 

-0.025(39) 

0.078(76) 

-0.007(153) 

-0.002(35) 

0.067(70) 

0.054(37) 

-0.033(38) 

0.051(37) 

-0.006(54) 

0.063(105) 

-0.067(38) 

-0.041(40) 

0.068(58) 

-0.015(29) 



-0.254(389) 

-0.173(30) 

-0.022(89) 

0.038(110) 

-0.090(18) 

-0.121(49) 

-0.104(21) 

-0.064(21) 

-0.056(20) 

-0.058(26) 

-0.182(79) 

-0.095(21) 

-0.026(18) 

-0.068(28) 

-0.038(13) 



0.149(14) 

0.126(39) 

0.150(8) 

0.134(11) 

0.071(23) 

0.131(9) 

0.119(11) 

0.132(9) 

0.119(13) 

0.115(9) 

0.116(17) 

0.143(29) 

0.096(14) 
0.106(17) 
0.076(13) 



-0.019(243) 

-0.007(29) 

0.038(43) 

-0.110(94) 

0.000(24) 

0.007(40) 

0.066(25) 

0.005(28) 

0.038(23) 

-0.046(44) 

0.070(63) 

0.019(16) 

-0.018(29) 

0.051(33) 

0.007(22) 



TABLE XV: Results for the isosinglet unpolarized generalized form factors for dataset 4. 



-°10 



AU + d 
^20 



yu+d 



/^u + d 



A'd+d 
■^30 



r}U+d 
-"30 



-t[GeV' 



A 



0.000 
0.177 
0.238 
0.240 
0.378 
0.459 
0.464 
0.498 
0.665 
0.738 
0.858 
0.875 
0.996 
1.042 
1.173 
1.235 



3.000(5) 
1.231(560) 
2.030(48) 
1.767(225) 

1.503(55) 

1.328(202) 

1.325(89) 

1.157(85) 

1.067(101) 

0.948(119) 



0.522(23) 

■1.770(2.373) 0.355(168) 



-0.140(211) 
-0.592(516) 

-0.013(184) 

-0.943(508) 

0.049(201) 

0.040(208) 

0.140(212) 

-0.136(266) 



0.426(18) 
0.378(43) 

0.367(18) 
0.360(52) 
0.331(22) 
0.298(30) 
0.311(28) 
0.244(37) 



-0.324(575) 

-0.026(62) 

0.049(125) 

0.013(61) 

0.144(122) 

-0.036(57) 

0.103(73) 

0.090(64) 

-0.014(100) 



0.153(25) - 

-0.461(825) 0.189(112) 0.033(574) 



-0.210(58) 0.154(16) 
-0.133(136) 0.145(20) 



-0.097(33) 

-0.090(93) 

-0.080(40) 

-0.074(38) 

0.003(39) 

-0.081(49) 



0.865(86) -0.173(175) 0.315(36) -0.048(66) -0.093(38) 



0.686(113) 
0.593(65) 



-0.033(235) 
0.073(118) 



0.196(35) 
0.212(23) 



-0.073(70) 
-0.065(42) 



-0.041(39) 
-0.043(21) 



0.129(19) 
0.115(22) 
0.098(16) 
0.070(25) 
0.111(21) 
0.073(36) 



0.070(21) 
0.092(22) 



-0.028(55) 
0.055(80) 

-0.082(45) 

0.012(73) 

-0.054(40) 

-0.087(56) 

-0.013(43) 

-0.066(83) 

0.057(31) 

-0.005(46) 
-0.017(37) 



TABLE XVI: Results for the isosinglet unpolarized generalized form factors for dataset 5. 



35 



^20 



20 



4 u+d 



-t[GeV^ 



0.000 
0.107 

0.124 
0.125 
0.219 
0.244 
0.245 
0.254 
0.359 
0.379 
0.471 
0.473 
0.508 
0.578 
0.615 
0.632 



3.000(6) 

2.610(460) 

2.437(35) 

2.465(229) 

2.445(482) 

2.067(42) 

2.113(213) 

1.956(91) 

1.763(54) 

1.644(91) 

1.541(72) 

1.601(234) 

1.468(81) 

1.318(68) 

1.306(104) 

1.257(73) 



3.374(1.535) 

-0.087(285) 

-0.247(499) 

1.221(1.054) 

-0.102(238) 

-0.505(455) 

-0.355(327) 

-0.237(229) 

-0.475(274) 

-0.274(261) 

0.283(552) 

-0.035(253) 

-0.286(204) 

-0.386(290) 

-0.104(188) 



0.510(24 
0.516(89 
0.468(20 
0.458(40 
0.497(93 
0.412(19 
0.422(40^ 
0.405(24 
0.366(20 
0.363(24 
0.357(24 
0.372(52 
0.362(30 
0.324(24 
0.292(25 
0.315(24 



-0.379(418) 

-0.071(90) 

0.071(139) 

0.025(268) 

-0.045(71) 

-0.112(124) 

-0.058(89) 

-0.061(69) 

-0.069(77) 

-0.120(79) 

0.097(148) 

-0.117(83) 

-0.060(65) 

-0.057(79) 

-0.062(60) 



-0.766(673) 

-0.334(76) 

-0.156(202) 

-0.560(221) 

■0.119(41) 

-0.250(103) 

-0.061(62) 

-0.091(44) 

-0.077(48) 

-0.117(49) 

-0.151(80) 

-0.101(61) 

-0.100(32) 

-0.112(43) 

-0.055(33) 



0.164(26) 
0.139(50) 
0.131(18) 
0.140(20) 
0.159(41) 
0.110(19) 
0.130(20) 
0.137(22) 
0.094(21) 
0.101(17) 
0.102(23) 
0.120(24) 

0.100(22) 
0.086(18) 
0.093(24) 



-0.325(359) 

0.091(71) 

-0.060(90) 

-0.094(165) 

0.043(52) 

0.025(84) 

0.072(67) 

0.066(52) 

0.055(58) 

0.021(71) 

-0.011(91) 

0.101(42) 

0.001(53) 

-0.022(52) 

0.066(55) 



TABLE XVII: Results for the isosinglet unpolarized generalized form factors for dataset 6. 



-t[GcV^] 



"^3(1 



r>u— 
^30 



0.000 
0.203 
0.243 
0.243 
0.420 
0.475 
0.477 
0.498 
0.697 
0.741 
0.911 
0.918 
0.996 
1.117 
1.199 
1.239 



1.141(12) 
0.963(91) 
0.973(12) 
0.930(37) 

0.786(80) 
0.849(13) 
0.780(37) 
0.831(18) 
0.742(17) 
0.713(18) 
0.683(24) 
0.595(45) 
0.658(17) 
0.575(23) 
0.575(27) 
0.584(17) 



9.251(403) 

9.159(853) 

4.361(1.326) 

6.514(227) 

5.515(570) 

6.798(305) 

4.682(232) 

4.604(222) 

3.913(286) 

3.175(512) 

2.690(176) 
2.806(225) 
2.866(147) 



0.280(4) 

0.244(22) 

0.256(4) 

0.234(10) 

0.216(21) 

0.231(4) 

0.208(10) 

0.226(5) 

0.211(6) 

0.203(6) 

0.203(8) 

0.173(15) 

0.198(6) 

0.176(8) 

0.171(9) 

0.181(6) 



0.740(136) 

0.562(286) 

0.303(343) 

0.497(71) 

0.390(161) 

0.539(107) 

0.432(81) 

0.445(81) 

0.486(116) 

0.241(204) 



0.100(2) 

0.087(10) 

0.093(2 

0.086(4 

0.081(9 

0.085(2 

0.078(4 

0.084(2 

0.079(3 

0.077(3 

0.081(4 

0.070(6 



0.254(67) 0.070(3 
0.366(93) 0.070(4 
0.216(68) 0.068(3 



0.066(71) 

0.090(120) 

0.129(186) 

0.063(39) 

0.017(78) 

0.188(55) 

0.094(46) 

0.164(48) 

0.117(64) 

0.006(100) 

0.052(37) 
0.149(51) 
0.068(43) 



TABLE XVIII: Results for the isovector polarized generalized form factors for dataset 1. 



-t[GeV^] 



0.000 
0.200 
0.242 
0.243 
0.414 
0.473 
0.475 
0.498 
0.693 
0.741 
0.904 
0.911 
0.996 
1.106 
1.195 
1.238 



1.136(18) 

1.115(124) 

0.967(16) 

0.974(51) 

0.902(107) 

0.848(16) 

0.819(49) 

0.828(23) 

0.736(22) 

0.710(24) 

0.691(32) 

0.738(85) 

0.621(22) 

0.580(28) 

0.565(37) 

0.559(20) 



9.891(469) 

10.370(1.088) 

6.917(1.552) 

6.654(253) 
6.279(672) 
6.481(360) 
4.549(262) 
4.508(266) 
3.778(377) 
4.214(829) 

2.733(228) 
2.342(288) 
2.640(168) 



0.278(6) 

0.236(27) 

0.253(5) 

0.237(12) 

0.217(26) 

0.231(6) 

0.220(13) 

0.224(7) 

0.214(8) 

0.199(7) 

0.191(11) 

0.170(24) 

0.189(7) 

0.173(10) 

0.174(12) 

0.173(7) 



20 



0.562(177) 

0.277(298) 

-0.199(414) 

0.442(96) 

0.391(195) 

0.383(127) 

0.380(90) 

0.278(106) 

0.232(169) 

-0.070(306) 



A u—d 

0.100(3) 

0.097(14) 

0.095(3) 

0.087(5) 

0.082(11) 

0.087(3) 

0.080(5) 

0.085(3) 

0.085(4) 

0.075(4) 

0.081(6) 

0.075(11) 



r>w — 
-°30 



-0.028(93) 

0.007(146) 

-0.091(246) 

-0.053(50) 

0.061(91) 

0.033(74) 

-0.015(67) 

0.025(60) 

-0.085(94) 

-0.022(161) 



0.212(92) 0.070(5) -0.031(54) 
0.110(130) 0.067(6) -0.004(76) 
0.187(85) 0.071(4) -0.076(55) 



TABLE XIX: Results for the isovector polarized generalized form factors for dataset 2. 



36 



-t[GeV'] 


Au—d 


0.000 


1.171(18) 


0.193 


0.926(117) 


0.241 


0.967(16) 


0.242 


0.934(49) 


0.403 


0.784(99) 


0.469 


0.854(18) 


0.472 


0.827(52) 


0.498 


0.811(23) 


0.685 


0.739(24) 


0.740 


0.711(26) 


0.891 


0.592(28) 


0.901 


0.669(79) 


0.996 


0.579(22) 


1.089 


0.562(31) 


1.189 


0.551(37) 


1.238 


0.529(22) 



r>u— 
-"10 



9.646(437) 

9.656(1.001) 

5.562(1.440) 

6.667(257) 

6.440(675) 

6.543(317) 

4.501(237) 

4.575(247) 

2.686(304) 

3.465(693) 

2.414(208) 
2.328(246) 
2.445(160) 



AU — d 

0.280(7) 
0.212(27) 

0.250(5) 

0.238(12) 

0.185(25) 

0.223(6) 

0.214(14) 

0.222(7) 

0.207(8) 

0.205(8) 

0.176(10) 

0.172(23) 

0.180(8) 

0.170(10) 

0.166(12) 

0.167(7) 



TJU — 

-"20 



0.601(166) 

0.684(315) 

0.581(359) 

0.484(73) 

0.495(189) 

0.479(109) 

0.446(75) 

0.505(92) 

0.205(142) 

0.315(242) 



0.101(3) 

0.068(15) 

0.092(3) 

0.089(5) 

0.074(11) 

0.083(3) 

0.079(6) 

0.085(3) 

0.078(4) 

0.079(4) 

0.061(6) 

0.060(11) 



0.102(89) 

0.171(139) 

0.324(245) 

0.067(49) 

0.044(101) 

0.094(68) 

0.068(68) 

0.100(62) 

-0.102(87) 

-0.068(129) 



0.317(75) 0.061(5) 0.039(47) 
0.297(107) 0.067(6) 0.117(65) 
0.297(72) 0.061(4) 0.033(51) 



TABLE XX: Results for the isovector polaxized generalized form factors for dataset 3. 



-t[GcV^ 



AU — d 



Su—d 
10 



AU — d 

^20 



riu — d 
^20 



AU 



Bu—d 
30 



0.000 
0.188 
0.240 
0.241 
0.394 
0.465 
0.469 
0.498 
0.678 
0.739 
0.880 
0.892 
0.996 
1.072 
1.183 
1.237 



1.162(32) 

0.931(204) 

0.945(26) 

0.934(80) 

0.710(166) 

0.815(27) 

0.894(86) 

0.797(36) 

0.740(37) 

0.696(38) 

0.564(51) 

0.898(206) 

0.574(33) 

0.471(45) 

0.519(71) 

0.475(28) 



9.410(668) 

9.110(1.700) 

4.255(2.373) 

6.498(379) 

7.418(1.076) 

6.362(470) 

4.611(376) 

4.407(354) 

2.530(492) 

4.339(1.694) 

1.935(280) 
2.172(415) 
1.999(229) 



0.257(10) 
0.224(49) 

0.228(8) 

0.222(21) 

0.184(43) 

0.214(9) 

0.218(23) 

0.195(10) 

0.210(11) 

0.192(12) 

0.156(16) 

0.232(56) 

0.174(12) 

0.144(15) 

0.139(22) 

0.152(11) 



0.357(248) 

-0.132(487) 

-0.140(627) 

0.457(135) 

0.616(293) 

0.297(184) 

0.495(133) 

0.570(152) 

-0.081(218) 

-0.070(658) 



0.091(6) 
0.082(31) 

0.088(5) 

0.078(9) 

0.087(21) 

0.086(5) 

0.078(10) 

0.078(5) 

0.082(7) 

0.072(6) 

0.062(12) 

0.041(21) 



0.034(152) 

-0.411(227) 

0.741(418) 

0.236(82) 

0.100(154) 

-0.006(102) 

0.217(115) 

0.245(102) 

0.011(151) 

-0.385(290) 



0.096(105) 0.058(8) 0.089(73) 
-0.011(177) 0.056(12) 0.023(118) 
0.131(120) 0.071(6) 0.038(80) 



TABLE XXI: Results for the isovector polarized generalized form factors for dataset 4. 
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37 



-t[GeV'1 



TJU — 

-°10 



TJU — 

-"20 



AU- 

^30 



30 



0.000 
0.177 

0.238 
0.240 
0.378 
0.459 
0.464 
0.498 
0.665 
0.738 
0.858 
0.875 
0.996 
1.042 
1.173 
1.235 



1.208(60) - 0.244(19) - 0.090(11) - 

0.557(398) - 0.172(96) - 0.189(90) - 

0.954(46) 11.827(1.139) 0.227(14) -0.112(427) 0.099(10) -0.019(302) 

1.105(152) 14.903(3.016) 0.239(37) 0.651(848) 0.088(19) 0.389(475) 



0.810(48) 
0.974(168) 

0.712(58) 
0.665(65) 
0.622(72) 
0.624(109) 

0.486(64) 



6.659(597) 
7.646(1.750) 

5.767(709) 
4.181(577) 
3.650(657) 
3.380(767) 



0.481(106) 2.296(589) 
0.408(50) 1.486(335) 



0.208(15) 
0.214(43) 

0.203(18) 
0.201(24) 
0.208(25) 
0.164(35) 

0.158(25) 



0.430(217) 
0.533(527) 

0.488(287) 
0.500(245) 
0.475(298) 
0.134(318) 



0.177(37) 0.305(285) 
0.137(18) 0.114(199) 



0.082(10) 
0.095(23) 

0.095(10) 
0.090(16) 
0.094(14) 
0.072(28) 



-0.087(168) 
0.539(307) 

0.301(187) 
0.397(246) 
0.201(212) 
0.346(249) 



0.067(21) 0.041(197) 
0.059(11) 0.090(147) 



TABLE XXII: Results for the isovector polarized generalized form factors for dataset 5. 



-t[GoV" AX 



Su—d 
10 



^20 



rtu—d 
^20 



A'" 



Bu—d 
30 



0.000 
0.107 
0.124 
0.125 
0.219 
0.244 
0.245 
0.254 
0.359 
0.379 
0.471 
0.473 
0.508 
0.578 
0.615 
0.632 



1.187(64) 

1.203(223) 

1.070(53) 

1.150(117) 

1.333(295) 

1.037(55) 

1.055(113) 

0.983(66) 

0.935(56) 

0.903(64) 

0.841(66) 

0.896(181) 

0.815(68) 

0.813(64) 

0.751(84) 

0.776(60) 



17.234(1.595) 

19.686(3.278) 

16.441(5.731) 

13.098(1.010) 

13.658(2.080) 

12.490(1.361) 

9.315(870) 

8.869(919) 

7.041(992) 

7.785(2.274) 

5.751(647) 
5.191(847) 
5.401(642) 



0.252(20 
0.262(51 
0.237(17 
0.251(28 
0.276(58 
0.230(16 
0.250(28 
0.218(19 
0.217(16 
0.214(19 
0.223(22 
0.260(53 
0.203(23 
0.204(20 
0.174(24 
0.192(19 



0.623(598) 

2.059(839) 

1.696(1.095) 

0.808(299) 

1.879(558) 

0.918(431) 

0.755(268) 

1.001(360) 

0.900(390) 

1.557(725) 



0.083(12) 
0.064(37) 

0.083(10) 
0.080(13) 
0.086(28) 
0.080(10) 
0.091(14) 
0.088(10) 
0.074(10) 
0.087(11) 
0.067(14) 
0.064(26) 



0.551(391) 

-0.454(561) 

-0.340(878) 

0.369(186) 

0.745(381) 

0.471(326) 

0.108(238) 

0.403(283) 

-0.100(289) 

-0.415(419) 
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TABLE XXIII: Results for the isovector polarized generalized form factors for dataset 6. 
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-t[GeV'1 



AU+d 

^10 



r>u+d 
-°10 



AU+d 
^20 



-°20 



Au+d 
^3(1 



r>u+d 
-°30 



0.000 

0.203 
0.243 
0.243 
0.420 
0.475 
0.477 
0.498 
0.697 
0.741 
0.911 
0.918 
0.996 
1.117 
1.199 
1.239 



0.623(17; 
0.497(79' 

0.537(14 
0.485(31 
0.425(75 
0.477(15 
0.414(38 
0.458(18 
0.413(19 
0.388(21 
0.404(25 
0.344(50 
0.365(19 
0.351(24 
0.327(31 
0.328(18 



5.234(429) 

4.941(958) 

2.380(1.559) 

4.086(272) 

3.397(726) 

3.924(375) 

2.776(283) 

2.681(305) 

2.561(305) 

2.056(605) 

1.945(214) 
1.522(279) 
1.711(189) 



0.170(6) 
0.138(22) 

0.158(5) 

0.144(10) 

0.136(22) 

0.145(6) 

0.122(11) 

0.142(6) 

0.134(7) 

0.127(7) 

0.125(9) 

0.110(18) 

0.124(7) 

0.113(9) 

0.114(11) 

0.114(7) 



0.657(161) 
0.676(308) 
0.624(461) 
0.545(107) 
0.304(243) 
0.559(144) 
0.430(111) 
0.469(128) 
0.429(127) 
0.282(240) 



0.061(3) 

0.057(14) 

0.058(3) 

0.052(4) 

0.056(10) 

0.055(3) 

0.047(5) 

0.054(3) 

0.054(3) 

0.049(3) 

0.048(5) 

0.042(8) 



0.297(92) 0.046(4) 
0.357(129) 0.044(5) 
0.308(107) 0.047(4) 



0.186(96) 

0.046(152) 

0.210(277) 

0.096(54) 

0.014(113) 

0.198(79) 

0.142(67) 

0.184(74) 

0.120(85) 

0.035(115) 

0.073(50) 
0.141(68) 
0.072(62) 



TABLE XXIV: Results for the isosinglet polarized generalized form factors for dataset 1. 



-t[GcV^ A 



u-\-d 



■DU-\-d 
-"10 



-"20 



A u-\-d 
^30 



B 



u-\-d 



0.000 
0.200 
0.242 
0.243 
0.414 
0.473 
0.475 
0.498 
0.693 
0.741 
0.904 
0.911 
0.996 
1.106 
1.195 
1.238 



0.627(23) 
0.663(107) 
0.531(18 
0.568(44 

0.543(94 
0.480(20 
0.535(47' 
0.468(24 
0.426(26 
0.404(25 
0.353(34 
0.383(86 
0.308(27; 
0.308(32' 
0.329(42 
0.279(26 



5.042(498) 

5.822(1.170) 

2.378(1.841) 

4.096(332) 

4.970(889) 

4.024(467) 

2.977(332) 

2.785(339) 

1.709(420) 

1.679(936) 

1.375(285) 
1.266(381) 
1.679(259) 



0.169(8) 
0.125(29) 

0.157(7) 

0.146(14) 

0.122(29) 

0.145(8) 

0.131(15) 

0.143(8) 

0.133(9) 

0.116(9) 

0.122(12) 

0.113(28) 

0.121(10) 

0.111(12) 

0.101(16) 

0.112(10) 



0.471(205) 

0.163(373) 

-0.367(566) 

0.392(129) 

0.255(294) 

0.315(183) 

0.298(136) 

0.050(159) 

0.278(183) 

-0.079(339) 



0.063(4) 

0.065(18) 

0.062(4) 

0.053(7) 

0.048(13) 

0.054(4) 

0.047(7) 

0.055(4) 

0.049(5) 

0.044(5) 

0.056(7) 

0.062(15) 



0.165(111) 0.050(7) 
0.036(180) 0.044(8) 
0.259(130) 0.050(5) 



0.073(123) 

0.053(179) 

0.096(365) 

0.008(66) 

0.012(124) 

0.108(102) 

-0.027(91) 

-0.043(87) 

0.040(115) 

0.172(177) 

0.135(64) 

0.060(102) 

0.057(78) 



TABLE XXV: Results for the isosinglet polarized generalized form factors for dataset 2. 
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-t[GeV^] 



^10 



r>u+d 
-°10 



^20 



ryu+d 
-°20 



AU+d 
^3(1 



rtu+d 
^30 



0.000 

0.193 

0.241 
0.242 
0.403 
0.469 
0.472 
0.498 
0.685 
0.740 
0.891 
0.901 
0.996 
1.089 
1.189 
1.238 



0.580(24) 
0.526(105) 

0.515(18) 
0.514(43) 
0.450(97) 
0.431(20) 
0.352(48) 
0.443(22) 
0.346(25) 
0.367(27) 
0.300(32) 
0.390(81) 
0.295(25) 
0.258(32) 
0.286(41) 
0.260(23) 



5.666(508) 

7.032(1.094) 

4.042(1.708) 

3.438(295) 

2.286(768) 

3.463(362) 

1.935(296) 

2.274(330) 

1.278(328) 

2.274(733) 

0.988(228) 
1.023(311) 
1.149(202) 



0.181(8) 

0.131(31) 

0.158(7) 

0.153(13) 

0.132(29) 

0.138(7) 

0.150(16) 

0.141(7) 

0.123(9) 

0.129(9) 

0.101(11) 

0.116(26) 

0.107(9) 

0.099(11) 

0.110(15) 

0.099(9) 



0.409(199) 

0.299(371) 

0.480(491) 

0.245(108) 

0.592(275) 

0.102(158) 

0.208(122) 

0.130(146) 

-0.112(154) 

0.087(271) 



0.072(4) 

0.046(20) 

0.062(4) 

0.064(6) 

0.041(15) 

0.055(4) 

0.063(8) 

0.062(4) 

0.051(5) 

0.056(5) 

0.036(8) 

0.042(14) 



0.071(114) 

0.017(183) 

-0.056(362) 

0.016(67) 

0.123(155) 

0.140(96) 

0.058(91) 

0.082(88) 

-0.255(123) 

-0.178(151) 



0.050(104) 0.037(7) -0.066(71) 
0.076(149) 0.049(8) 0.098(89) 
0.084(124) 0.038(5) 0.034(75) 



TABLE XXVI: Results for the isosinglet polarized generalized form factors for dataset 3. 



-t[GcV^ 



B 



u+d 



AU + , 
^20 



^u+d 



r>u+d 
-"30 



0.000 

0.188 
0.240 
0.241 
0.394 
0.465 
0.469 
0.498 
0.678 
0.739 
0.880 
0.892 
0.996 
1.072 
1.183 
1.237 



0.598(38) 

0.393(194) 

0.474(30) 

0.467(71) 

0.522(177) 

0.427(33) 

0.419(81) 

0.394(36) 

0.407(45) 

0.433(47) 

0.336(56) 

0.463(164) 

0.301(39) 

0.297(53) 

0.213(77) 

0.248(34) 



4.326(808) 

6.043(1.796) 

5.751(2.875) 

3.081(455) 

3.280(1.215) 

3.088(543) 

2.207(462) 

2.993(499) 

1.348(555) 

3.079(1.375) 

0.992(358) 
0.526(503) 
1.051(280) 



0.159(14 
0.177(56 

0.143(11 
0.134(22 
0.172(52 
0.136(12 
0.148(27 
0.124(12 
0.133(15 
0.126(15 
0.108(21 
0.178(57 
0.109(15 
0.097(20 
0.096(28 
0.083(14 



0.200(289) 

-0.240(548) 

0.146(821) 

0.076(174) 

0.494(432) 

0.316(264) 

-0.143(200) 

0.263(226) 

-0.253(281) 

0.320(604) 



0.064(8) 

0.088(40) 

0.059(7) 

0.054(11) 

0.049(26) 

0.057(7) 

0.064(12) 

0.052(7) 

0.045(8) 

0.049(8) 

0.051(15) 

0.043(29) 



-0.008(204) 

-0.291(284) 

-0.101(625) 

0.106(110) 

0.118(225) 

0.099(145) 

-0.026(158) 

0.135(142) 

-0.119(208) 

-0.369(374) 



-0.081(163) 0.035(11) 0.021(105) 
-0.094(262) 0.052(15) -0.030(158) 
0.126(174) 0.042(8) 0.051(112) 



TABLE XXVII: Results for the isosinglet polarized generalized form factors for dataset 4. 
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-t[GeV^1 



^10 



r>u+d 
^10 



4n + d 
^20 



^20 



AU+d 
^30 



-°30 



0.000 
0.177 

0.238 
0.240 
0.378 
0.459 
0.464 
0.498 
0.665 
0.738 
0.858 
0.875 
0.996 
1.042 
1.173 
1.235 



0.519(68) 

0.003(440) 

0.469(51) 



7.022(1.388) 



0.604(146) 8.708(3.368) 



0.321(56) 
0.536(170) 

0.340(62) 
0.202(75) 
0.128(85) 
0.312(124) 

0.234(74) 

0.337(124) 
0.172(54) 



2.498(733) 
3.063(2.143) 

2.307(873) 
1.395(771) 
-0.480(819) 
1.348(923) 



1.644(738) 
-0.034(422) 



0.123(26) 
-0.092(122) 
0.126(19) 
0.094(42) 

0.113(22) 
0.113(56) 

0.092(22) 
0.090(30) 
0.057(29) 
0.109(46) 

0.063(31) 

0.095(44) 
0.053(22) 



-0.146(558) 
0.199(1.052) 

0.069(314) 
0.183(767) 

0.114(422) 
0.180(357) 
-0.626(415) 
0.313(388) 



0.626(361) 
-0.405(295) 



0.049(16) 
0.268(131) 

0.065(14) 
0.031(25) 

0.053(13) 
0.033(29) 

0.056(13) 
0.056(19) 
0.052(18) 
0.114(41) 



0.072(29) 
0.046(16) 



0.182(426) 
-0.209(618) 

-0.143(228) 
-0.064(436) 

0.088(290) 
0.648(331) 
-0.066(277) 
0.618(338) 



0.168(254) 
0.163(218) 



TABLE XXVIII: Results for the isosinglet polarized generalized form factors for dataset 5. 



-t[GcV^ A 



u-\-d 



^10 



AU-'rd 
^20 



^20 



4^ 



u-\-d 



B. 



u-\-d 



0.000 
0.107 
0.124 
0.125 
0.219 
0.244 
0.245 
0.254 
0.359 
0.379 
0.471 
0.473 
0.508 
0.578 
0.615 
0.632 



0.484(72) 
0.732(209) 

0.493(54) 

0.722(109) 

0.867(246) 

0.500(58) 

0.544(118) 

0.503(67) 

0.478(63) 

0.427(73) 

0.379(74) 

0.558(178) 

0.456(75) 

0.399(76) 

0.348(93) 

0.433(62) 



11.813(1.845) 

16.409(3.651) 

10.112(5.404) 

7.939(1.175) 

6.702(2.649) 

8.214(1.553) 

6.080(1.065) 

4.748(1.289) 

3.233(1.004) 

5.559(2.237) 

2.998(772) 
2.366(973) 
4.459(837) 



0.143(25 
0.184(59 

0.139(20 
0.154(34 
0.226(61 
0.141(20 
0.176(38^ 
0.123(25 
0.147(22 
0.128(26 
0.122(25 
0.153(57 
0.110(30 
0.113(24' 
0.083(33^ 
0.113(26^ 



0.832(668) 

0.730(1.118) 

0.704(1.398) 

1.075(372) 

1.910(815) 

0.844(662) 

1.257(379) 

1.757(534) 

0.847(451) 

0.722(811) 



0.040(15) 
0.045(41) 

0.049(14) 
0.026(18) 
0.046(35) 
0.048(13) 
0.067(18) 
0.038(13) 
0.050(13) 
0.045(14) 
0.068(21) 
0.057(36) 



0.673(521) 

-0.717(732) 

0.104(1.269) 

0.348(271) 

1.500(526) 

0.195(499) 

0.540(343) 

0.294(434) 

0.784(383) 

0.284(551) 



0.725(283) 0.051(17) 0.462(205) 
0.515(464) 0.039(20) -0.024(295) 
0.926(457) 0.031(15) -0.426(344) 



TABLE XXIX: Results for the isosinglet polarized generalized form factors for dataset 6. 
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